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1 Setting the problem 



An nxn Toeplitz matrix Tn is a matrix with coefficients of the form {Tn)jk = Cj-k, < j, k < n — 1, 
for some given sequence {Ql^ez- An nxn Toeplitz determinant Dn(T) is the determinant of 
some given Toephtz matrix T = T^. li ip = (/^(e*^) is an integrable function on the unit circle 
= {z = e*^} (oriented in the positive direction) with Fourier coefficients 



ip[ e 



de 

2^' 



the Toephtz matrices Tn{ip) and Toephtz determinants Dn{^) associated with ip are given by 

(1) = {V'i-fc}o<j,fc<n-l 

and 

(2) DM = det r„((^) 

respectively. Note that if = c^, then T„ is self-adjoint. This is true, in particular, if is real 
valued, and it follows that for such functions Lp, Dn{y:>) and the eigenvalues A^"''', . . . A^"''' of Tn{^p) 
are real. 

Toeplitz matrices and determinants are named for Otto Toeplitz, who, in his Habilitationsschrift 
in 1907 (see |Toepl| |Toep2| ), initiated the study of quadratic forms Tiipj^XjUk with coefficients of 
special type (pjk = 'Pj-k- At the time Toeplitz was a young researcher in Gottingen where Hilbert 
was developing his general and abstract theory of functional analysis. Toeplitz introduced the forms 
T,Lpj_kXjyk, together with the associated matrices Tn{p) = {pj-k}Q<j^k<n-i and determinants 
Dn{'p) = det(T„(99)), in order to give concrete examples of Hilbert's general theory which could 
be analyzed in great detail. Toeplitz called the forms S Pj-k xj yk "L-forms" because of their 
connection, via Fourier theory, with Laurent series '^'^^o 'fk ■ 

A striking example of the connection between Laurent series, and more broadly analytic function 
theory, on the one hand, and L- forms on the other, is given by the following problem studied by 
Caratheodory: For which sequences cq, ci, . . . , of complex numbers does the polynomial 



p{z) = Cq + CiZ-\ h CnZ" 

have an extension to an analytic function p{z) in the unit disk ID 

ip{z) =do + diz-\ h (i„ + dn+i 2;"+^ 



(3) 



= {\z\<l} 
+ ... , 



dj = Cj, < j < n 



such that ^p{z) > in B? In |Carl| in 1907, Caratheodory characterized such sequences in terms 
of Minkowski's theory of convex bodies, but a few years later in 1911 Toeplitz |Toep3| gave the 
following algebraic characterization (see also |Car2j . |CarFej| ): p{z) = cq + ci z + ■ ■ ■ + Cn z^ has an 
extension p{z) as above with ^p{z) > in B if and only if the Toeplitz matrix 
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with c_fc = Cfc, < k < n, is non-negative definite. Over the years, up till the present, there 
have been many applications of this beautiful result, both theoretical and applied (see e.g. | Weg| ) . 
Many of the results and applications in this paper, however, rely on the connection between Toeplitz 
matrices/determinants and another problem in analytic function theory, the Riemann-Hilbert Prob- 
lem (RHP) — see, in particular, [CG| |LSj and the references therein. We note that Wiener-Hopf 
factorization theory and the Wiener-Hopf method fall under the general rubric of RHP. 

An extraordinary variety of problems in mathematics, physics and engineering, can be expressed 
in terms of Toeplitz matrices and determinants. This is a remarkable fact that could hardly have 
been anticipated in 1907, and there is a vast literature on the subject. A classical reference is the text 
by Grenander and Szegqj [GreSz] in which the authors describe applications of Toeplitz matrices to 
problems in analytic function theory, stationary stochastic processes and linear estimation theory 
in statistics. In addition to problems in probability theory and statistics. Fisher and Hartwig 
[FisHartT] present a long and varied list of applications of Toeplitz determinants to problems 
in statistical mechanics, including the ground-breaking work of Kaufman and Onsager on the 
Ising model. Many results in the theory of Toeplitz determinants can be found, for example, in 
the seminal papers of Harold Widom (see |Widl] |Wid2] |Wid3] |Wid4| . amongst many others) 
and in the books of Bottcher, Silbermann, and Grudsky (see |BottSilbl| [BottSilb^ [BottSilbS] 
|BottGr| , amongst others) . For a review of some more recent developments in the theory of Toeplitz 
determinants, see |Krj . 

Closely related to Toeplitz matrices are Toeplitz operators, defined as follows (see e.g. [Doug] , 
[BottSilbS] ). Let denote the unit circle as above with Lebesgue measure 27 ; and let L'^{S^) 
denote the associated Hilbert space. The Hardy space is defined as the closed subspace 

|/ G L\S') : ^ f{^'') de = Q, n = 1, 2, . . .| . 

The orthogonal projection Pj^ of L'^{S^) onto is given by 

00 00 

-00 
Now let 93 S L°°{S^). Then the Toeplitz operator T{ip) : — )• associated with 99 is given by 
(5) T{^)f = P+vf, feH\ 

The operator T{ip) is bounded in H^, and with respect to the standard basis {e*'^^ : A; > O}, T{ip) is 
represented by the semi-infinite matrix acting in i'^ = {u = (uq, ui, . . . ) : Yl'o' < 

00}. Thus the Toeplitz matrices Tn{ip) are finite sections of the Toeplitz operator T{ip). The 
function ip is frequently referred to as the symbol oi the operator T{lp) and the sequences {T„((/?)}„>i, 
{D„((^)}„>1. 

Matrices with coefficients of the form {cjj^k}^^. ^<;„_^ are called Hankel matrices. These ma- 
trices are in turn finite sections of the semi-infinite Hankel operators with matrix coefficients 
{^j+k}j fc>o acting in The theory of Hankel matrices and operators is parallel to, and closely 
related to the theory of Toeplitz matrices and operators (see e.g. (jl66p et seq. below). Operators 

^For many years "Szego" was mis-spelled in the literature as "Szego" . Knuth first introduced the Hungarian umlaut 
into TeX specifically so that names like Szego, Erdos and Konig would be spelled correctly. (Private communication 
to the authors from D. Knuth.) 
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acting in L^(0, oo) with kernels of the form K(x — y), or of the form K{x + y), x,y > 0, may be 
viewed as continuum hmits of Toephtz matrices, or Hankel matrices, and are called Wiener-Hopf 
and Hankel operators respectively (see e.g. [BottSilbS] |Pelj ) . 

Our goal in this paper is not to review the above applications and developments in a systematic 
way, but rather to focus on one central problem in the theory of Toeplitz determinants, viz., 
the Szego Strong Limit Theorem (SSLT) (see Theorem 3 below). We will, in particular, trace the 
history of this theorem and its later generalizations, in response to a succession of questions raised 
in the analysis of the Ising model. An earlier, and very informative, description of the influence of 
the Ising model on SSLT and related developments, was given by Bottcher in [Bottlj . 

The story begins in 1915 when Szego [Szlj . at the very start of his career, proved the following 
result conjectured earlier by Polya. 

Theorem 1. Let (/^(e*^) > be a continuous, positive function on the unit circle and let Dn{(p) 
be the associated Toeplitz determinant. Then 

(6) hm - log Dniif) = ^ r logv(e'') dO. 

Note that Dn{f) > 0. Indeed for u = [uq, mi, . . . , n^)^ S C"^"*^, a simple calculation shows that 

•^"'^ fc=0 ^ 

and as ^{e^^) > 0, we conclude that Tn{ip) is strictly positive definite, T„(</?) > 0. Thus Dn{^p) = 

det Tn{(p) > 0. In addition, the eigenvalues A^"'^ . . . , aI"'^ of Tn{ip) are positive and so ([6]) can be 
rewritten in the form 

.ogAS" + ... + logAi" 1 r,^^ ^ 

Szego recognized that d?]) is a special case of equidistribution for the A^"^ in the sense of Weyl, and 
in 1920 [Sz2j Szego proved the following result (see also jGreSzj ). 

Theorem 2. Let Lp[e^^^ be a real-valued function in L°°{S^) with m < (/^(e'^) < M a.e. If F{\) is 
any continuous function defined on the interval m < X < M , we have 



Taking F{\) = log A, we clearly recover d?]). We consider more detailed asymptotic properties 
of the eigenvalues below (see Section [9|). 

Relation ([6]) may be rewritten in the form 



(9) L'„(v3) = exp 



n j log 99(^6*'') 1^ + 



as n — )• 00. The issue of determining the precise nature of the term o(n) came to the fore in the 
following way. 
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The two-dimensional Ising model, the central model in statistical mechanics (see [McWul] |Palj 
|Lie| |Domlj ). named for E. Ising (see [Is]), concerns the interaction of random spins aij = ±1 at 
sites in Z^. Of great interest is the situation in which only nearest-neighbor spins interact 

and the interaction energy is given by 

— Jl CTjj CTjj + l — J2 CTjJ CTj+lj 

where the vertical and horizontal interaction constants, Ji and J2 respectively, are translation 
invariant, and in addition, Ji and J2 are positive, so that the system is ferromagnetic, i.e. parallel 
spins have lower energy than anti-parallel spins. As is standard in statistical mechanics, one analyzes 
the infinite system on by first considering the spins in a finite rectangular box A C Z^ of size 
M X N (see below for a discussion of boundary conditions) and then letting M, N ^ 00. For such 
a box the total interaction energy of the spins a = {(Jij) is given by 

(10) E\{a) = - ^ [Jl aij CTjj+i -|- J2 fJjj crj+ij] 

A 

where Ji cjjj- Cij+i is included in the sum if or (i, j -|- 1) G A and J2 cjjj- Cj+ij- is included if 
(i, j) or (i -|- £ A. The associated normalized Gihhs measure is then 

(11) P^^(^) = ±e-^^('^)/'=^^ 

where is Boltzmann's constant and T is the temperature. Here the partition function Zp^ is 



(12) Zk = Y1 



where the sum is taken over all possible spin configurations a in A. For fixed finite sets A C A, one 
defines the correlation functions 

(13) ( n =e n ^m-p^a(^)- 

The thermodynamic limits of these correlation functions as A f Z^, i.e. M, N ^ 00, 
(14) 




are the objects of principal physical interest. In particular such correlations \^n(j j)eA ^i-i/ should 
model physical phenomena such as phase transitions. For example, a bar magnet has a critical 
temperature Tc called the Curie point: for temperatures T < Tc, the magnet exhibits spontaneous 
magnetization, but for T > Tc the magnetization is zero (in zero external field). This behavior 

cannot be described by a function such as ^0(1 j)eA (fT3|) . which is real analytic for T > 0: 

in the thermodynamic limit, however, analyticity can be, and indeed is, destroyed. 

In the presence of an external magnetic field, one must add a term of the form h j^^^ aij to 
Ej^{a), /i G M, but the most complete results have been obtained only when the field is absent, and 
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we will always assume h = 0. It remains an outstanding problem in the theory of the Ising model 
to find the analogs of the explicit formulae below for the correlation functions in the case h ^ 0. 

We note the following: In order to compute E\{a) in (jlOp . we need to know cjj.j not only at the 
points G A, but also at the points adjacent to A. In particular, if M and are integers and 
A = G : < i < M, < j < AT}, then we need to know ctm+i,j, < j < N and 

ai^N+i, Cj^-i, <i < M. In practice, this issue is addressed by imposing boundary conditions on 
A. One common choice is, for example, periodic boundary conditions, ai^M+i,j = ^nj, crij+N+i = 
(Jij. Another choice is the fully magnetized boundary conditions where cjjj = +1 (alternatively 
cTij- = —1) at all points adjacent to A. General boundary conditions are discussed, for example, 
in [Gal| . and also from a different point of view, in [LebM-Lofj . For each choice of boundary 
conditions, provided the thermodynamic limits in (jl4p exist, we obtain a family of functions 



In the spirit of the classical moment problem, one says that the family in (jlSh determines a (Gibbs) 
equilibrium state (see e.g. [Siml] [Ruj ) . A priori, different boundary conditions can give rise to 
different families (jlSp . and hence to different equilibrium states. For example (cf [Pal] ) . for closed 
boundary conditions aij = +1 adjacent to A, one finds (cto,o) > 0, but for periodic boundary 
conditions, by symmetry, (cto,o) = 0. Nevertheless, it suffices for our purposes to note that for the 
two-point correlation function that we discuss below, the thermodynamic limit is independent of 
the boundary conditions (cf |BGJ-LS ]). The same is true for the free energy F per unit spin in 
(fTTl) below. 

In 1 dimension, Ising [Is] showed that the Ising model does not exhibit a phase transition for 
any temperature T > 0. What about 2 or 3 dimensions? In a landmark paper in 1936, Peierls [Pel] 
gave a simple argument asserting that in 2 or 3 dimensions, the Ising model does indeed exhibit 
spontaneous magnetization at some temperature Tc > 0. It turns out that Peierls' argument 
involved an incorrect stepH, which was corrected only many years later in 1964 by Griffiths [Grif] . 
Griffiths confirmed Peierls' conclusion and "Peierls argument" remains a standard tool in many 
statistical situations. The first exact quantitative result for the 2-dimensional Ising model was 
obtained by Kramers and Wannier in 1941 [KraWan] when they wrote down the following formula 
for Tc. In the case Ji = J2 = J, 



In 1942, Wannier gave a talk on this work at a meeting of the New York Academy of Sciences. In 
an extraordinary denouement, at the end of the talk Onsager announced that he had obtained an 
exact solution for the two-dimensional Ising model (without magnetic field). Onsager published 
the proof of his result two years later in [Onslj . More precisely, Onsager had obtained a formula 
for the partition function Za, and, in particular, in the case of a square lattice A^r of size N x N, 
with Ji = J2 = J, he showed that in the thermodynamic limit the free energy per unit spin F 

•^First discovered by M. Fisher (see [Grif] ). 



(15) 





(16) 
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exists and is given by 



1 



F = -ksT lim — log Zaj^ 



(17) =-kBT 
= -kBT 



log ( 2cosh ^— + ^ y #1 / (i</'2log (1 - -(C0S(/>1 + COS(/>2; 







where k = 2 sinh (2 J/Zc^r) / cosh^ (2 J/fc^T). It follows from this formula that the specific heat 
C = —T diverges logarithmically as T —t- Tc it 0, where Tc is precisely the critical temperature 
obtained by Kramers and Wannier (note that for Tc we have k = 1, and the argument of the 
logarithm in the double integral vanishes at the integration limit (pi = 4>2 = Onsager's analysis 
allows for different interaction constants Ji / J2, in which case equation (fT6]l for Tc must be 
replaced by 

(18) sinh ( sinh ( = 1. 



ksTcJ \kBTc 

For the description of the events that followed, we refer to the informative, and entertaining, 
recent papers by Baxter |Bax| and McCoy |McCl| . In the mid-1940's, Kaufman, a doctoral student 
at Columbia University, began working with Onsager, who was then at Yale, on the 2-D Ising 
model. Onsager's computation of the partition function had stunned the physics community by its 
ingenuity and complexity, and Kaufman was determined to find a simpler approach. As Onsager 
relates |0ns2) , Kaufman "... decided to look for a possible connection with spinor theory. Why 
not? In fact, it seemed like very good sense, and so it was .... By the summer of 1946 she had a 
beautifully compact computation of the partition function, bypassing all tedious details." 

The methods of Onsager and Kaufman for the partition function are based on the calculation 
of the eigenvalues of the transfer matrix [KraWan] associated with the problem. Another simplified 
solution based on this approach was given in 1964 by Schultz, Mattis, and Lieb |ScML] . On the 
other hand, it was observed in 1941 by van der Waerden [vanderW] that the calculation of the 
Ising partition function is equivalent to the problem of counting closed polygons on the square 
lattice. This observation was the origin of an alternative (and more direct) combinatorial approach 
to the Ising model. A solution using such an approach was given by Kac and Ward in 1952 
[KacWj (see also [PotWar] where the calculation of the correlation functions we discuss below is 
addressed by this method). Kac and Ward did not intend to present a rigorous derivation; the 
missing proofs were given later by Sherman [Sherj and Burgoyne [Burj following a conjecture by 
Feynman. A rigorous combinatorial solution by a related method was given by Hurst and Green 
|HurGr| . Another combinatorial solution was obtained by Kasteleyn [Kasl] . A useful modification 
of Kasteleyn's approach is due to Fisher |Fis3j (see [McWul] ) . One more simple solution based 
on the method of Kac and Ward was given by Vdovichenko [Vdov| (see iLanLifj ) although [Vd QVj 
omits the discussion of the boundary conditions. In [Kasl], Kasteleyn maps the Ising problem onto 
the dimer covering problem (i.e. the problem of counting all possible coverings by dimers) on a 
special lattice, slightly more complex than the square lattice. It is interesting to note that the dimer 
problem on the square lattice, which thus can be regarded as a simpler version of the Ising problem, 
was solved by Kasteleyn himself [Kas2| and, independently, by Temperley and Fisher |TemFis) in 
1961 (see also [LieLoss] for a later alternative, particularly simple, solution of the dimer problem 
on the square lattice). 
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We now return to 1946. Kaufman's method and results were eventually published three years 
later in [Kaul] . Onsager continues in |Qns2| : "By itself that was only a more elegant derivation of 
an old result; but the approach looked powerful enough to produce a few new ones. Very well, how 
about correlations?" As noted by Baxter, the history of the Ising model from that time forth has 
been the study of these correlations. 

In 1948, at a conference at Cornell, and again in 1949 at a conference in Florence, Onsager 
astounded the audience by announcing that he and Kaufman had recently obtained an exact formula 
for the spontaneous magnetization Mq of the 2-D Ising model. He gave the result as 

(19) Mo = (1 - C.) , k^ns - (sinh ^ sinh g,) ' 

for < /cons < 1, corresponding to T < Tc. And for fcons > Ij corresponding to T > Tc, Mq = 0. 

Onsager and Kaufman never published a proof of this result. The situation was as follows. 
In 1949 they published their famous paper [KauOns] . based on Kaufman's spinor (free-fermion) 
approach, in which they found, in particular, an explicit expression (see (j2ip below) for the two- 
point correlation function along db row in the tlicrmociyriaiiiic limit, i.e. {o'\^i(7i^i-\-n 

). By physical 

arguments 

(20) Mo = lim 

n— ^-oo 

SO that in order to derive ()19p . what they had to do was to control the limit in (|20p . As Onsager 
describes the situation [Ons3j : "And then, finally, it became apparent that we had a last problem, 
of the degree of the order, and that turned out to need more invention." 



2 Toeplitz meets Ising 



What precisely was this "last problem" that Kaufman and Onsager had to face in order to compute 
the limit in (I20p ? It is here that we make contact with the theory of Toeplitz matrices. In [KauOns] 
Kaufman and Onsager expressed (cri^i cri.i+„,) as a sum of two Toeplitz determinants as follows (in 
the notation of [MPW| ): 



(21) 



where 
(22) 
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bo 
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■ bi-n 


1)" = c^2^et 
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bo 
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■ b2-n 
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bn-3 
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. b.J 




( bl 


b2 


b3 .. 
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■ bn-1 
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C2 = cosh K2, S2 



sinhi^^, 



tanh 



J2 
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(23) '^^^l Jo "^"^^ ^ '^^^^ ~ ^ 



(24) cos 6{9) = sin {6* (d)) sin 6 cosh 2K^ - cos {6* (d)) cos d 

and v'ons(e*^) = e*^*^^^ is a function introduced previously by Onsager in |Onsl) 



(25) V'onsfe 



i9 



1 — 7i e*^ \ / 1 — 72 e 



1 — 7i e *^ / \ 1 — 72 e 



where 71 = 2:3, 72 = 2^2/^1 

(26) ,1.-^1 . "^2 . I-Z2 

z\ = tanh - — — , zo = tanh - — — , Zn = , 

ksT' ksT' 2 ^^^^^ 

The branch in (j25|) is chosen so that V3ons(e*'^) = e*'^*^'^-' > 0. 

Remark 2.1. As noted in [MPW] . there is a sign error in formula (j2ip in [KauOns] . 

So the problem that Kaufman and Onsager had to face to compute Mq via ()20p . was the purely 
mathematical problem of computing the asymptotics of n x n Toeplitz determinants as n — )• 00. At 
the time, however, all that was known was Szego's result ([9]) with unknown error term o(n). But 
to compute (I20p . this error term is precisely what one needs to know! 

Note that 

(27) < 71 < 1 and 71 < 72. 
If 72 7^ 1 , we see that 

(l-7ie^^)(l-72e-^) 



(28) V^ons e 



1(1-71 e-*«) (1-72 6 



and hence (/7ons(e*^) is a smooth, single- valued function on the unit circle S^. On the other hand, 
if 72 = 1, then 

(29) V^ons (e^')=^ e"^^/^ j^—^ , < ^ < 2vr. 



1 1 - 71 e 

Simple algebra shows that for Xi = Ji/ksT, i = 1,2, 

72 ^ 1 <^ sinh [xi + X2) - cosh (xi - X2) ^ 

and 

sinh^ (Xi + X2) - cosh^ (xi - X2) = k',^^ - 1 
where fcons = (sinh 2xi sinh 2x2) ""^ as in ([T9|). Thus 
(30) 72 ^ 1 ^ A:ons ^ 1. 
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In other words, sub-critical (resp. super-critical) temperatures, T < T^. (resp. T > T^) correspond 
to Z2 < z\ (resp. z\ > z\). And, of course, T = Z2 = zi. 

It turns out that in their initial (unpubhshed) calculations to compute Mq, Kaufman and 
Onsager did not actually use ([2T]). What happened is described by Onsager in |0ns2| and |0ns3] 
(see also [Bax] for more details, particularly concerning the Wiener-Hopf calculation below). In 
addition to (f2T]) . Kaufman and Onsager had derived, but apparently did not publish (see private 
communication to C. Domb |Dom2l p. 201]), an expression for the 2-point correlation function 
along a diagonal (cJi.i, o"i-|_„^ which is much simpler than (I2ip . Only one Toeplitz determinant 
is involved, 

(31) {cr 1,1 a i^n, l+n) = Dn (V'diag) 

wherq3 



1 



(32) yJdiag 



iQ\ I ^ "'ons C \ 



e 



"'ons c 

and again A;ons is given as in ()19p and ¥?diag(e*'^) > 0. The same physical arguments as in ()20p yield 

1 

(33) Mo = lim (fii 1 (Ti+„,i+ri) ^ . 

n— ^00 

Onsager then realized that the eigenvalue equation for the Toeplitz matrix r((y9diag), 

^ n— 1 

lifc = 1" ^ (V5diag)fc_j- lij, < A; < n - 1 

was a discrete analog of the Milne integral equation 

(34) f{x) = A{x- y) f{y)dy, x > 



which arises, with a particular, explicit function A{-), in radiative equilibrium theory (see, for 
example [DymMc] ) . Milne's equation can be solved using the Wiener-Hopf method, a technique 
with which Onsager was familiar from lectures that Wiener had given earlier in the Mathematics 
Department at Yale. So Onsager tried the Wiener-Hopf technique on the above eigenvalue equation 
for T((/3diag), to obtain the eigenvalues A = A^"'', I < £ < n, and then computed Dn{ipdiag) = 
Y[e=i ^e^^- Letting n — )• 00, formula (|19p for Mq then emerged. This was the basis for the first 
announcements of the result in 1948 and 1949. So why didn't Kaufman and Onsager publish the 
details of their calculation? The situation was as follows. As in the case of Milne's equation, the 
Wiener-Hopf method requires very detailed knowledge of the kernel {^Pdiag)j-kj in the case of Ising. 
Onsager sensed that there was another prize to be had: What about Toeplitz determinants with 
general symbols 93? And so, before Kaufman and Onsager could get around to publishing their 
calculation, Onsager began looking for a method to evaluate asymptotically Toeplitz determinants 
with general symbols. We quote Onsager [ Ons3] : " ...and lo and behold I found it. It was a 
general formula for the evaluation of Toeplitz matrices. The only thing I did not know was how 



^In |Dom2| . the symbol is (pdiag, not <^diag- However, the determinant is the same, -Dn((^diag) = detTn(v?dia 

detr„((pdiag)^ = detT„((pdiag)^ = detT„(((3diag)^ = detT„(((9diag) = -D„ ((ySdiag) • 
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to fill out the holes in the mathematics and show the epsilons and deltas and all of that, and the 
limiting processes; I did not know just how it should be done and what mathematicians really knew 
about limiting processes in that ball park." As it turned out, by that point, mathematicians knew 
a lot! Before it was clear what conditions to place on the symbol, Kaufman and Onsager spoke to 
Kakutani and then Kakutani spoke to Szego. Szego then revisited his calculations from 1915, and 
in 1952 [Sz3j evaluated the o(n) term in ([9]) for a very general class of symbols (p. This celebrated 
result is the Szego Strong Limit Theorem mentioned above. Faced with Szego's general result, 
Kaufman and Onsager published neither their first, nor their second method. As Onsager noted 
|0ns2| . " . . .the mathematicians got there first." In the long history of mathematics and physics, 
it is most unusual for a physicist to be scooped out of a formula by a mathematician! 

There is, however, more to the story. On the web (see |Baxj ) there is the draft of a paper 
from 1950 (never published), almost certainly by Onsager and Kaufman, that describes their more 
general method and obtains the result (jl9p for Mq. Here the authors use formulae ()20p and ()2ip . 
There is also an earlier letter [Kau2j written on April 12, 1950, from Onsager to Kaufman containing 
some of the calculations from the above draft. Presumably to allay any concerns about priority, 
Onsager assures Kaufman at the end of the letter that "There will be time for all these things." 
As one reads [_0ns2J |0ns3] , one senses the faint hint of regret that Onsager must have felt for such 
optimism. In a letter from Kaufman to Onsager, written on May 12, 1950, a month after the letter 
|Kau2] ■ Kaufman says "Here is a draft of Crystal Statistics IV". In |Baxj . Baxter makes a strong 
case that the above paper from 1950, now posted on the web, is indeed Kaufman's draft of Crystal 
Statistics IV. 



ike 



log(^(e*^) dO, A; G Z, denote the Fourier coefficients o/ log (/?(e*^) . Then 



3 Szego's Theorem 

Here is Szego's result |Sz3j [GreSzj 

Theorem 3 (Szego Strong Limit Theorem). Let 99(6*^) be a positive, C^^^ , e > 0, function on . 
Let {\ogip)k = ^jl^e 

(35) 
where 

(36) 



E{^) = J2k\ilogip)f. 



k=l 



□ 



Thus the o(n) term in d 
Remark 3.1. Note that 

(37) 



is given by E{ip) + o(l). 



g(logv)o =gxp 



— / log^(e'"]de 



is the geometric mean of ip. 

In contrast to (j35]) . Szego's earlier result Theorem [Tl limn-^c>o {Dn{ip))" = e(i°g'^)o, is known 
simply as Szego's First Theorem, or sometimes just Szego's Theorem. 
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In the years that followed, there was considerable effort by many mathematicians to weaken 
the assumptions in Theorem [3j Along the way, many new methods, quite different one from the 
other, were discovered to prove ([35|) . These developments are described in the outstanding recent 
monographs by Barry Simon [Sim2| and Bottcher and Silbermann [BottSilbS] : also see [McC2] 
and (Bottlj for earlier summaries. In |Sim2j . the contributions of Kac [Kac] . G. Baxter [BaxGlj . 
Hirschman |Hir| . and Devinatz [Dev| . in particular, are discussed, leading up through successively 
weaker assumptions to the definitive result of Ibragimov [Ibrj in 1968 giving necessary and sufficient 
conditions on (p[e^^) for (I35p to hold under the assumption of the positivity of ip. 

Theorem 4 (Ibragimov). Let (/?(e*^) ^ be a probability measure on and suppose that log(/7(e*^) 
is integrable. Then (jSSp is always true in the following sense: lim„_i.oo gi^ioit^ ^.^''J^o.ys exists and 
equals e^^'^\ including the case where one, and hence both, are infinite. □ 

One can consider Toeplitz matrices Tn{dfi) for general probability measures dfi{6) = (/?(e*^) ^ + 
dfj,s{0) on S^, where dfis denotes the singular part of dfi. We have Tn{dfi) = {fJ-j-k}o<j,k<n-i, where 
fi£ = f e~^^^ dfj,{9), ^ G Z and = det Tn{dfi). The following remarkable result generalizes 

Theorem [TJ 

Theorem 5 ( [SzT] [SiS] [Sil] ) . Let dfi{e) = (^(e*^) §^ + djJsiO) be a probability measure as 
above. Then Szegd's Limit Theorem is always true in the following sense: lim„_j>oo Dn{dfi)" always 



exists and equals exp f^^log(p[e^^) ^ 
□ 



, including the case where one, and hence both, are zero. 



Szego proved the result when dfis = 0, whereas Verblunsky was able to handle the case dfis / 0. 
Of course, as (^(e*^) S L^ (f^); XT^ log ^ < and the integral can diverge, possibly, only 

to — oo. The striking feature of Theorem [5] is that lim„^oo ™ is independent of d^s- 

Is it possible that (f35]l remains true when d/is 7^ 0? Here the definitive result is due to Golinskii 
and Ibragimov. 

Theorem 6 ( |GolIbrj ). Let dfi{6) = (p[e^^) ^ + d^s{0) be a probability measure as above, and 
suppose logip{e'^) G L^ {^) . If dfis / 0, then 



(38) lim 



exp 



n/_^log(^(e^^)f 



+00 



□ 



In other words, if ([55]) were to hold with E{ip) = X^fcLi ^ K^OS v)fc|^ < co, then necessarily 
dus = 0. 

In [Sim2j . Simon presents six different proofs of (j35p . Some of these proofs require stronger 
conditions on ip (see [Sim2| ). We now describe these methods in some detail (but without giving 
precise conditions on ip in each case) in order to illustrate the extraordinary variety of mathematical 
areas and techniques that inter-relate with Szego's Strong Limit Theorem. Our presentation follows 
|Sim2j . After this presentation, we will also describe three additional proofs not covered in [ Sim2) . 
together with some brief comments on Szego's original proof (see Remark 13. 2p . 

The first proof uses ideas and results from the theory of orthogonal polynomials on the unit 
circle (OPUC's). These polynomials were introduced by Szego in the early 1920's [Sz4| |Sz5| in the 
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course of his investigation of the eigenvalues . . . ,X^^ of Toephtz matrices T„. The classical 
references for OPUC's are |Sz5| [Gej . The OPUC's Pk{z) = Xfc + . . . , Xfc > 0, associated with a 
probability measure dfi{9) on are formed by orthonormalizing 1, z, , . . . , , . . . with respect 
to d^l{e) in L2 (51^ diJL{d)), 



(39) / ^ Vk{e'')V3 (e'') di^iO) = 6kj, k,j>0. 



s 

OPUC's are intimately related to Toeplitz determinants. For example 

where = — (Xj+i)^^ Pj+i(0)) j ^ Oi are the so-called Verblunsky coefficients. The ^j's 
have modulus less than 1, and so F{dfi) = liiRn^oo Dn+i{dfi) / Dn{dfi) always exists and equals 

Xo^rij^o ~ l^il^)- -^^^ then lim„^oo {Dn{df^))" always exists, and has the same limit. So the 
proof of Szego's Limit Theorem, boils down to showing that 



(41) A n (1 - l^.f) = exp r log^(e^^) d9 



where diJ,{9) = <f{e''^) ^ + d^s{9). In turn, one finds that 
(42) lim^^limn--^^-'^^'^)"" 



n~.^{F{d^,)r — nr=„-i(i-ie 



and the proof of SSLT for d^s = reduces to showing that the limit on the RHS is just e^^"^-*. 
As we will see below, OPUC's play a crucial role in analyzing the Fisher-Hartwig conjecture for 
Toeplitz determinants with singular symbols (^(e*^). 

The second proof utilizes a remarkable identity which gives an exact formula for Dnif), 

(43) DM = e"(i°g'^)o+^('^) det (1 - Q„ H{b) H(S) Q„) , n>l 

where det denotes the determinant in £2 = £2{'^+), Z+ = {0, 1, . . . }, Qn is the orthogonal projection 
in £2 onto = {u = (uq, ui, . . .) ^ £2 : Uj = for < i < n}, and H(b), H(c) are the Hankel 
operators with kernels j>o, {c-i-j-i}ij>o- Here 6(e*^) = ^jbje''^^, c(e*'') = ^jCje"-^^ 

are defined in terms of the Szego function V{z) = exp J^^ fw^f log (^(e*^) d9 



(44) 6.^, 

where V = I?(e*^) = lim^_^gie |^|^]^ ^(z). The point is that H{b)H(c) is a trace-class operator in 
£2 = ^^0' hence Qn H{b) H(c) Qn — )• in trace-norm by abstract arguments as n — 00. But 
then det (1 — Qn H{b) H(c) Qn) — )• 1 as n — )• 00, by the continuity properties of the determinant. 
(For more information about the trace class and determinants, see [BottSilbS] [Sim S].) Thus SSLT 
is immediate, once one has proved (143 p . Formula ()43p is known as the Borodin-Okounkov formula 
and was proved [BorOk] in 2000 in the context of combinatorics and random matrix theory. The 
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formula evoked great interest, and meanwhile several new proofs and generalizations have been 
given (see, e.g., [BasWidj [BottWidT] and |Sim2j ). However, it turns out that it was already 
proven many years earlier by Geronimo and Case in 1979 in their work on inverse scattering theory 
[GerCase] , The broader significance of this paper was not appreciated at the time. But more to the 
point, in a separate section in their paper titled "Szego's Theorem", Geronimo and Case actually 
used (H3]l to prove SSLT. It is unfortunate that these developments were overlooked by the experts 
in the field. 

The third and fourth methods utilize the following Heine-type multi-integral representation for 
£>„((/?) (see, e.g., |BottSilb3) [Sim2]) 

(45) DM =^,r---r n - ^"H' n ^ (^"0 f • 







0<i<A:<n-l j=0 



In the third method, due to Bump and Diaconis [BumpDia] , the authors observe that the RHS 
of (j45p can be re-written via Weyl's integration formula as Jm^^^ e^^ia) dg^ where dg denotes Haar 

measures on the unitary group U{n), and Fn{g) = X]j=o lo§ '/'(^*^^*'^'') > where je^^-'^^-'} are the 
eigenvalues of g. Substituting the Fourier expansion log(/9(e*^) = Yl'k=-Qo^^^'^)k e*'^^, one obtains 



/ exp (^(log(/.)fctr(/) I 

Ju{n) J 



(46) L»„(^) = e"(^°s^)° 
and after expanding out the exponential, one ends up with the representation 

(47) L»„((^) = e"('°s^)° 



Here the sums are over all fe-tuples of non-negative integers t = {ti, . . . , t^), s = (si, . . . , s^), and for 
any Tjts, rjs's are exphcit constants depending on {(log 99)^}. 

The magic of the method is that one can use the theory of characters together with Frobenius- 
Schur-Weyl duality to compute Tt{g) Ts{g) dg explicitly, obtaining 



(48) / Ttig)Tsig) = 6tsW{t) 

Ju{n) 

where W{t) = nj=i i*^ ^i'' provided Yl^=i3^j Sj=i 3 — ^- Substituting (jiSj) into ([T7|) 

and controlling the errors as n — )• c«, one immediately obtains (j35p . 

In the fourth proof, due to Johansson [Johl| . one rewrites (j45p with 92(6*^) = e^^^^ in the form 

(49) Dje^) = A L-^H„ie,,...fi„-^) JK-i 
^ ' ^ ' n\ J 27r 27r 

where 

n-l 

j=0 0<A:<j<n-l 
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and 

(51) 



w{e) = - log 



This expression exhibits n!Z)n(e^) as the partition function of a gas of 1-dimensional particles 
with 2-dimensional Coulomb interactions — a "log-gas" in the terminology of Dyson. Dyson 
introduced this terminology in the context of his analysis |Dyl| of Circular Ensembles of n x 
n random matrices. For such ensembles Dyson found that the distributions of the eigenvalues 
|g«ej |n-i^ < 0j < 27r, of the matrices were given by 



(52) 

Here /3 
(53) 



1 



■'n,l3 



n 



' 27r 



da 



n-l 



0<j<k<n-l 

1 for the orthogonal, (3 = 2 for the unitary, and /3 

2^ " 



2tt 



4 for the symplectic ensemble, and 



Z. 



n,l3 











7n-l 

2tt 



n 

0<j<k<n-l 



is the normalization constant. Dyson then noted that the distribution of the eigenvalues {e*^-'}?=o 



was identical with the distribution of the positions < < 2tt, j = 0, 



1 of charges 



in a finite gas with 2-D Coulomb interactions at inverse temperature (5 and with Gibbs mea- 



sure ^ eM-PHn{eo, On-i)}^ 



2lT 



where Hn{Oo 



7n-lj 



E 



0<j<k<n- 



_ilog|e*' 



6*^*^1, and the partition function is given by (j53p . As emphasized by Dyson [Dyl| , a consequence 
of this identification is that ". . .the thermodynamic notions of entropy, specific heat, etc., can be 
transferred from the Coulomb gas to the eigenvalue series. This will prove very useful, as it gives 
us a precise and well-understood language in which to describe the statistical properties of the 
eigenvalues" . If the matrices U in the ensembles carry a weig ht e-^*'''3('^), then the energy in the 
corresponding Coulomb gas must be replaced by 



n-l 



(54) 



Hr, 



E 

0<j<fe<n-l 



log|e" 



j=0 



where Q now plays the role of an external force acting on the charges. 

In operational terms, the identification implies, in particular, (see |Dy 2|) that the free energy 
F = —j3~^\ogZn,p is approximated to leading order as n ^ oo by the classical thermodynamic 
recipe, F ~ Ft, where 

(55) /3Ft = G2 + Gi + Go. 
Here G2 is the macroscopic Coulomb energy 

(56) G2 = -^j J a{e'^)a{e'^')log\e''^ -e'^'ldOdd' + ^ 
and Gi + Go is the local free energy term 



ie 



)de 



(57) 
(58) 



Gi 



Gn 



1 



a{e'%\oga{e'')de 



n 



/3\ , 27r , ^ 
^ log— + logr 
2 / n 



logT 
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The function a{e^^) > is the macroscopic density of the gas and is chosen to minimize the 
functional on the right-hand side of (j55p subject to the constraint 



(59) J a{e'^)d9 = n. 

Dyson first apphed these thermodynamic notions to compute |Dy2| the probabihty Pg that there 
are no eigenvalues for a random matrix in the interval (0, 2s /n), in the bulk scaling limit where the 
average distance between the eigenvalues is rescaled to 1, for the Gaussian Unitary Ensemble (see 
[Mehj ) ■ Using these notions, he derived, for the first time, the leading order of P^, ~ exp{— 
as s — )• oo. (We return to this problem further on (see (j22ip et seq).) A year later in [Dy3| , Dyson 
showed how these notions play out in the context of SSLT. First Dyson noted, as explained above, 
that the Heine-type representation (j45p immediately displayed n!D„(e^) as the partition function 
of a log-gas as in ()49p - (|5ip above. This connection between Toeplitz determinants and the Coulomb 
gas was first pointed out in an unpublished manuscript by G. Kreisel, which was made available to 
Dyson by courtesy of E. WignerO 

For (3 = 2 the term Gi = and the constant Gq = — logn!, and we have 

(60) Ft = -^J J a{e'^)a{e'^' )log\e''^ - e'^'\ded9' + J Q{e'^)a{e'^)de - ^lognl 
The constrained minimum is obtained when 

(61) Q(e^^) - [ (T(e*^')log|e^^ - e'^'\de' = const 



on the unit circle. Solving this equation together with (j59p . Dyson finds that for Q{e^^) sufficiently 
smooth, 

^ oo 

(62) FT = --\ogn\ + nQo- J2 \^\QkQ-k 

k=~oo 

where {Qk} are the Fourier coefficients of Q{e^^). For V = —(3Q = —2Q, we obtain 

1 1 °° 

(63) --log{nlDn{e^)) = F^ Ft = -- log n\ + nQo- ^ \k\QkQ-k 

k=—oo 

which is precisely SSLT, D„(e^) ~ expjnVo + ^jV^p}. There is an irony in this argument: 
on the one hand, SSLT was proved in response to a problem in statistical mechanics, but on the 
other hand, SSLT is obtained here by the methods of statistical mechanics. The above calculation 
is based on the (unproven) validity of the thermodynamic recipe (I55p ~ (j58p . (]59p : Johansson's proof 
in |Johl| in effect makes this approach rigorous. More precisely, the proof in [Johlj starts with the 
identity -Dn(l) = 1 and develops a perturbation theory around ip{z) = 1. Johansson notices that 
the minimum of the double sum in (I50p , and thus the maximum of the integrand in (I49p for V = Q 
is reached when the points e*^^, j = 0, . . . ,n — 1, are uniformly distributed on the unit circle, i.e. 
when they are the vertices of a regular n-gon. For ^{z) = e^^^-* ^ 1, normalize f by requiring 



*We note that already in the 1880's Stieltjes gave an electrostatic interpretation of the zeros of certain classical 
orthogonal polynomials (see [Sz5j . Section 6.7). For a discussion of logarithmic potentials from the viewpoint of 
potential theory, see the book of Saff and Totik [SaTo] . 
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Vq = 0. If one adds the corresponding first sum with V ^ to (jSOp . the position of the maximum 
of the integrand in (j49p becomes displaced, but in a sense only slightly if n is large. For large n 
Johansson introduces the following change of variables in 



(64) 



ipj h(e 

n 



i = 0, . . . ,n - 1, 



<1) 

logZ)„ (1 - A/i) = - V — tr {[Tn{h)r) 

^— ' m. 



m=l 



where h{z) = -iJ2'k=iiVkz''-V^ j^z ) is the conjugate function to V{z). The main contribution to 
the integral comes from a small neighborhood of a configuration where ijjj 's are uniformly distributed 
on the unit circle. Due to this crucial fact, all the sums which appear (as a result of the change of 
variables (jM]) ) in the exponent, apart from the main one, X^o<A;<j<n-i ^(V'fc ~ V'j)) can be replaced 
by integrals with a good precision and produce, as n — t- oo, the constant (j36p . while the remaining 

h I exp(- Eo<k<,<n-l Wii^k - i^,)) ^ = ^n(l) = 1. 

The fifth method is combinatorial and is due basically to Kac |Kacj . At the start one renormal- 
izes </3(e*^) so that supg ip[e'^^) = 1. Setting h = 1 — ip and w\ = 1 — Xh, < A < 1, we have wi = ip. 
Under the additional assumption that inig (/^(e*^) > 0, we clearly have < h{6) < supg h{6) < 1. 
Using the identity log det A = tr log A, one expands log Z)„ (1 — Xh) for < A < 1 in a power series 
(note that ||T„(/i)|| < 

(65) 

and also 
(66) 

On the other hand 
(67) 

and so 
(68) 

where 



log (1 - Xh) 



2^ 



m=l 



27r 



(log(l-A/i)), 



m=l 



k=l 



Y,k\{iog{i-xh))f = Y,s^ — 

fit 

m=. 



oo m— 1 



m 



k=i e=i ^ ^ 
Thus the strategy is to show 



'ike 



e'l'lfr-'iO)—) 
27r / 



(69) 



lim ■ 

n— >oo 



tr ([r„(/i)])"-n / h^ie) 



de 

2^ 



Combinatorial issues arise when one tries to control the complicated multilinear sums in tr [T„(/i)]™. 

The sixth method is based on a formula of Baik, Deift, McLaughlin, and Zhou (unpublished) 
for relative determinants 



(70) 



log 



1 

dt 
Jo 



2lT 



d 
It 



log ipt 



Kt{e) M9) 



de 

2^ 
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where (/? and ip are two functions on and ift = {l — t)+ tip and Kt{9) is the Christoffel-Darboux 
kernel 



n-l 



(71) Kt{9) = Y,\Pj{e'',Vti^ 

j=0 

and pj[e''^, (ftip) are the OPUC's associated with the weight 99t(e'^) ^/^(e*^) SSLT then follows 
by computing the (weak) limit of ^ Kt{9) as n ^ cc. Thus, as in the first method, SSLT is seen 
here as following from properties of OPUC's pj(e*^). 

The mathematical challenge that is posed by SSLT can be viewed as follows. In functional 
analysis one commonly considers continuous maps F, say, from one fixed topological space X to 
another fixed topological space Y. Then in order to verify that lim„,_j>oo F{xn) = F{x) all one has 
to do is verify that — )• x as n — )• oo in the topology of X. In particular, if X is the space of trace 
class operators on a separable Hilbert space %, and F{x) = det (1 + x), then F : X — )■ y = C and 
lim„_^oo det (1 + x„) = det (1 + x) if x„ — )• x in trace norm. But SSLT is not of this standard type. 
Instead, one has F : M(n, C) — )• C taking G M(n, C) to F{An) = det A„. So the initial space 
X = Xfi = M{n, C) is changing with n and it is not at all clear how to topologize the situation so 
that ^^liuin F {An) = F(lim^„)". If, for example, we let X = M(oo,C) be the inductive limit of the 
M(n, C)'s, then Tn{(p) converges to the Toeplitz operator T(ip) in the associated topology, but T{ip) 
is not of the form 1 + trace class, and so F(T{ip)) = det(T((/?)) is not defined. In such situations 
when X = Xn, or y = y„, is varying, there is no general procedure to follow and each problem 
must be addressed on an ad hoc basis. And when there is no preferred path to the solution of a 
problem, there are many paths, as we see from the six very different methods above. In this regard, 
the method of Borodin-Okounkov-Geronimo-Case is singled out: after factoring out e"(^°S'^)o"'"-^('^\ 
one has 

(72) -^^P^^^ = detil-Bn) 

where i?„ = Qn H{b) H(c) Qn is trace class in £^ and Bn — )• in trace norm as n — )• oo. In other 
words, one of the paths is to renormalize the problem so that it assumes the standard form. 

This pathway, to reduce SSLT to standard form, goes back to Harold Widom, who was the first 
to apply operator-theoretic methods to the problem of Toeplitz asymptotics. In |Wid4j . Widom 
gave yet another proof of (|35p . This proof, which was found more than 35 years ago, has not lost its 
charm. We follow the presentation in [Bottlj . It turns out that by a simple trick, the determinants 
of finite sections of the inverse Toeplitz operator T((/?)~^ are easy to compute, 

(73) det Pn T{ip)-^ Pn = e-"(i°g 

where Pn projects u = {uq,ui, . . . ,Un-i,Un, ■ ■ ■) G to (uq, ui, . . . , n„_i) G C". Simple algebra 
shows that 

detTni^p) = det PnT{ip)Pn = det(P„r-i {^-^) Pn) 

^^^^ X det(l + (P„r-1 (99-I) Pn)'^ PnKPn) 



where K = T{ip) — T ^ [tp is trace-class. By ([73]) . 

(75) detP„T-i((/j-i)P„ = e"('°*^^)o. 
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But (P„ T ^ {^{p ^) Pn) Pn can be shown to converge strongly as n — > oo in £^ to (T ^ ((/? ^)) ^ = 

T (v?""^), and as K is trace class, Rn = [PnT~^ ^n) ^ Pn K Pn ^ T {^'^) K. We conclude 

that as ra — )• oo 

(76) = + ^ + ^ (^"') ^) = ^ ^(^)- 

Using the remarkable Helton-Howe-Pincus formula (see [HeltHow] ) . det^e"^ e~^) = e*''[^'^l, 

if A, B are bounded, and [A, B] = AB — BA is trace class, one then shows that det(T (y?"^) T{ip)^ = 
e^(^), and 1^ is proved. 

Of all the methods to prove SSLT, perhaps the simplest and most direct was introduced by 
Basor and Helton in [BasHelt j . As in |Wid4j . the method also proceeds by reducing SSLT to 
standard form. We follow |BasHeltj : Suppose (p has a Wiener-Hopf factorization cp = (/?_|_93_, 
where (p+, (/?^^ are bounded analytic functions in {|2| < 1}, and ip-, ipZ} are bounded analytic 
functions in {\z\ > 1}. Then by direct calculation, the Toeplitz operator T{(p) has a factorization 
T{ip) = T{{p^)T(ip^), furthermore, 

(77) PnT{ip+) = PnTi^+)Pn, T{^^)Pn = P„T(v?_)P„. 

As T{(p±)T{(p^^) = 1, we see that T{(p±)~^ exist and T{ip±)'^ = T{(p^^). We thus have 

(78) PnT{ip)Pn = PnT{^+)PnT{ip+)-'T{^^)T{ip+)T{ip^)-^PnT{^^)Pn 

SO that 

Dniv) =det(P„r((^+)P„)det(P„{r(vp+)-i,T((^„)}P„)det(P„r(99_)P„) 
= ((9^+)o(^-)o)"det(P„{T(^+)-\r(9._)}P„) 

where \^T{ip^)~^ ,T{ip^)^ = T(ipj^)^^T{ip^)T(ipj^)T{ip^)^'^ is the multiplicative commutator. How- 
ever {r((/?+)^^, T(99_)} = 1 + T((/?+)^-'^ [T(99_), T((y9+)] T(99_)~^ and under mild smoothness condi- 
tions on ipie"^^), the commutator [T{(p^),T{ip^)] = T((/?_)T((/?+) — T{ip^)T{ip^) is trace class. For 
such (p, det(P„ {T((/9+)~^, r((y9_)| P„) is equal to the Fredholm determinant 

det(i + p„t(<^+)-i [r(<^_), r(v9+)] r(v9_)-ip„) 

which converges by abstract theory as n — )• oo to 

det(i + r((^+)-i [r((^_),r((^+)]r((^_)-i) =det{r(99+)-\r((^_)} 

(80) =det(r(v9+)-ir((^_)r(<^+)r((^„)-i) 

= detT{ip-^)T{ip). 

Again by the Helton-Howe-Pincus formula, we have 

(81) lim ^"^'^i = detT{p,-^)T{ip) = e^(^) 

which proves SSLT under mild assumptions on ip (see jBasHelt] Theorem 2.1.). 

In [Deilj . the author gives yet another proof of (135p by expressing Dnif) as a Fredholm deter- 
minant of the following form: 

(82) DM=detil-Vn) 
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where is a trace class operator acting on L'^[S'^) 

(83) Vn f{z) = I Vn [z, Z') f {z') dz' , \z\ = 1, 

J\z'\=l 



with 



The operator Vn is oscillatory and does not converge as n ^ oo in the trace-norm topology. 
However, Vn has the special form of an integrable operator, that is, operators with kernels of the 
form X^jLi "^^ ^ ' with z, z' on some contour F in C, m < oo, and fj, gj are given functions 
on r. Such operators were first singled out as a distinguished class by Its, Izergin, Korepin and 
Slavno^l (see [Deilj for a pedagogic presentation). Associated with the integrable operators there 
is a canonical Riemann-Hilbert Problem, which, in this case, is of oscillatory type. Such problems 
can be analyzed asymptotically by applying the nonlinear steepest descent method for Riemann- 
Hilbert Problems introduced by Deift and Zhou in |DZ] . and further developed in [DVZj . and SSLT 
then follows (see [Deilj ). 

Remark 3.2. In the methods of Widom and Basor-Helton the existence of the limit of ^i^iog'^^)^ as 
n ^ oo follows directly from general theorems in functional analysis: the identification of the limit 
with e^^'^^ is a separate matter and uses other means (Helton-Howe-Pincus) . The situation is similar 
regarding Szego's original proof of SSLT (see |Sz3j |GreSzj ): as ^p{e'^^) > 0, the ratio ^i^iog'^)^ is 
monotone increasing in n, so the existence of the limit (which may, a priori, be infinite) is immediate. 
Again the evaluation of the limit is a separate matter. Indeed, Szego first computes the limit for 
a dense class of symbols {^} for which Dn{(p) can be evaluated explicitly for n sufficiently large, 
n > n^p, and then obtains the general result by an approximation argument if ^ ip. 

By contrast, all the other methods we have described above prove the existence of the limit, 
and compute its value, both at the same time. 

We return now briefiy to the 2-dimensional Ising model. One notes that Szego's Theorem [3] 
does not actually apply to V'diag(^)) say, in (j32]) . because the symbol is not real valued. Moreover, 
Szego's proof of (f35]l depends in a crucial way on the positivity of 93 (e*^) (see Remark 13.21 above), 
and it was not at all clear at the time how to extend Theorem[3]to non-real symbols (p. So Onsager's 
lament that ". . . the mathematicians got there first" , was in fact not true! The result in the 1950 
unpublished draft of Kaufman and Onsager mentioned above, without "the epsilons and deltas and 
all of that", but allowing for non-real ip, was in fact the first to break the tape! Now although 
the symbol <^diag(^) is not real, for < /cons < 1 however, log (/?diag(^) is still smooth on S^, in 
particular, the winding number of </'diag(^) on is zero. The first result for such symbols was 
proved by G. Baxter in 1963 (BaxG2] . He showed that under certain technical conditions, if ¥?(e'^) 
is smooth, non-zero, and has zero winding number on S^, then 

(85) hm = e^^- Miogv^),{iog^)_,_ 

Note that for v'diag(^') ™ (|32|) . one has for < ko^s < (log (/^diag)^ = (sgn i) koisf^l^l for ^ 7^ and 

(log (/?diag)o = 0- Hence Y.eLiH^oS^)i (log (/?)_£ = -3 = 3 log (l - k"^^^), which yields 

(fT9D via (El. 



Some of the elements of the theory of integrable operators were already present in the earlier work [Sakh] . 
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In retrospect, many of the proofs of SSLT do extend directly to complex- valued symbols with 
smooth logarithm. This is true, in particular, for Kac's proof of (j35p in 1954 [Kacj, as noted by 
Montroll, Potts and Ward in [MPWj (see footnote 17, p. 316). Moreover, the proofs in |Wid4j and in 
the text |BottSilb3] . for example, are all designed from the very beginning to cover complex- valued 
functions. 

If (^(e*^) is non-zero and continuous, say, on S^, and has zero winding number, then f{e^^) = 
e^(<^' ) for some continuous, periodic function V{e^^) = log 99(6*^). The strongest version of SSLT 
with complex symbols ^p{e^^) is the following. 

Theorem 7. Let y(e*^) G L^{S^) be a (possibly complex-valued) function on with Fourier 
coefficients {Vkjj^^x satisfying 

00 

(86) 1^1 l^fcl'<°o ' 

fc=— 00 

then 

(87) = kVuV-. _ 

□ 

Note that for V{e^^) satisfying ([86]), e^^^""'^ G LP(5\ for 1 < p < 00 (see [Sml2] ). The proof 
of Theorem [7] is given by Johansson [Johl] . In particular, his argument for the extension of the 
result to complex-valued V{z) proceeds by considering the analytic functions 



and applying Vitali's Theorem. When z is real, \\m.n->oo gn{z) is given by (i85]l . and to obtain (|87|l 
one just sets z = i. 

Under the additional assumption that V is in L°°{S^), Theorem [7] was already established in 
|Wid4j . We also note that condition (j86p can be replaced by certain "non-symmetric" requirements 
on the decay of the Fourier coefficients Vk and V-k-, A; > 0; see, for example. Theorem 10.32 and 
Corollary 10.42 of |BottSilb3j . 

Remark 3.3. In |Johl| . Johansson observes the following very interesting consequence of Theorem 
[7] (see also [ Joh2| . |DJ) ). Let V{e'^^) be a real-valued function on with Fourier coefficients {Vk} 
satisfying (|86l) and suppose Vb = / V{e^^) If = 0. As in the Bump-Diaconis method for SSLT 
described above, equip the unitary group U{n) with Haar measure dg and consider the random 
variable Fn{g) = TJ^jZlV {e''^'^^'^) where (e^^^^^)} are the eigenvalues of G U{n). Then we have 
the following Central Limit Theorem: 

(00 
2^ k\Vk[ 
k=l 

where A^(0, cr^) denotes the normal distribution with mean zero and variance a. Indeed, for any 
real t we have 

(89) [ e''^"^3) dg = Dn{e''^) 

Ju{n) 
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and so by Theorem [71 as n — oo 

Ju{n) 

from which we conclude (j88p . One notes that as opposed to the standard Central Limit Theorem 
for independent random variables, the 0j's are very strongly uniformly distributed and no 
scaling is needed. 



4 Back to Ising 

After Onsager announced formula (jl9p for Mq in 1948, and again in 1949, but did not provide a 
proof, three years went by until a derivation of the formula was found by Yang in a tour de force 
in 1952 [Yanl] , Yang's method is based on results of Kaufman and Onsager, but his approach 
is different, and he does not explicitly use Toeplitz determinants. In his Selected Papers in 2005 
( [Yan2] ). Yang described the long and winding road that led him to the Kaufman-Onsager formula 
pgp . The story has much charm: 

"I was thus led to a long calculation, the longest in my career. Full of local, tactical tricks, 
the calculation proceeded by twists and turns. There were many obstructions. But always, after 
a few days, a new trick was somehow found that pointed to a new path. The trouble was that I 
soon felt I was in a maze and was not sure whether in fact, after so many turns, I was anywhere 
nearer the goal than when I began. This kind of strategic overview was very depressing, and several 
times I almost gave up. But each time something drew me back, usually a new tactical trick that 
brightened the scene, even though only locally. 

Finally, after about six months of work off and on, all the pieces suddenly fitted together, 
producing miraculous cancellation, and I was staring at the amazingly simple final result." 
Remark 4.1. Let £'a (<t) denote the interaction energy in a rectangular box A of size M x as in (jlOp 
above, and let E\^h{a) = E\{a) — ^ j)eA denote the energy in the presence of a magnetic 
field of strength h. Let Zx,h = Y^^e~^^'^^'^'>f^^'^ denote the associated partition function, and let 
Fh^h = ~if^ log.^A,h denote the associated free energy in the box A. What Yang proved in |Yanlj 
was that, for T < Tc, 

My = — lim lim ^'^/^ ^ 

^g-j^^j HO \A\~^oo h/N 



1 - ^L)^ =Mo= lim cJi,„)t 



1 

oo 



(In [Yanl J . Yang only considered the case Ji = J2- Yang's method was extended to the general 
case Ji / J2 by Chang [Cha] .) 

As noted and emphasized by Schultz, Mattis and Lieb [ScMLj in 1964, neither My nor Mq 
coincide, a priori, with the physically fundamental definition of the spontaneous magnetization 
Mphys, viz.. 



^phys = - lim lim ^ff'^ 
^ ^ HO |A|->oo oh 



(92) ^ 

= lim lim 



hiO |A|^oo MN 



j2 er'^ ^^'^^"^^ 
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Only some years later was it shown that indeed for T < T^. 
(93) Mphys = My = Mo = (1 - kl^,) 



(see, in particular, |LebM-Lo^ . |BGJ-LS] ). 

In 1955 an important technical advance was made by Potts and Ward |PotWar] when they 
showed, in particular, that the correlation function along a row for the 2-D Ising model 

could be expressed in terms of a single Toeplitz determinant, rather than a sum of two Toeplitz 
determinants as in ()2ip . They were unable at the time to verify that their solution agreed with 
(I2ip : This was done later in 1963 [MPWj by Montroll, Potts and Ward, who also used recent ideas 
of Kasteleyn on the Pfaffian approach to the dimer and Ising problems to give a new and simpler 
proof of the formula in { PotWar] . The formula in [PotWar] |MPW] is the following 

(94) (o-l,l = Dn{(Pons) 

where ¥?ons(e*^) is the Onsager function defined in ()25p . Through discussions with Mark Kac, 

Montroll et al. were aware of SSLT. Furthermore they knew, as noted above, how to extend Szego's 

result via Kac's method to the case where the symbol is complex. For T < Tc, Montroll, Potts and 

Ward noted that < 71 , 72 < 1 by ([30|) , and so (p"^^ 1 (e^^) / and 

ifons has no winding on S^. As 

iag( 



in the computation for (^diag(e*^) above, we find 



(95) (log<^ons), = -(^^-|j, i>l 

(96) (logV^ons)_,= (|-|), i>l 



and (logv3ons)o = 0. This implies Y."^! H^og (fons) e (logV3ons)_£ = I log J and after 



(i-7?)(i-7l) 
(1-7172) 

some elementary algebra using ([26|) . we again obtain p^ . but now via (pO|) (cf. the discussion 
following dSSD). 

We now consider the spontaneous magnetization for the Ising model for temperatures T > Tc. 
As fcons t 1 when T f T^, one anticipates from (jl9p that at T = the spontaneous magnetization 
Mo is zero. In [KauOns] . Kaufman and Onsager proceed in the following way. They derive an 
approximate formula for in (j23p (6^ = 'Y_i^ in the notation of [KauOns] ) at T = 

(97) ^'=~^r(2FTT) ^ 

They then take (j97p as an equality for all A; > 0, and substitute this relation into (j2ip to obtain an 
approximate asymptotic formula for (ui^i fii^i+n) in the form of a linear combination c^^ Wn — S2 Wn 
of two Cauchy determinants Wn and Wn = W-n- Such determinants can be evaluated explicitly in 
terms of gamma function^ 

(98) |ir„| = rn rb + i)rfa + i) 

'jr A. A. 



IT 



ir(g + i)r(g + i)' 



n-l 



(99) ^'^=3^nf 



2 'ff T{q + l)T{q + l) 



,ir(9-|)r(g + |) 



Formulae (|98p (I99[) differ slightly from the corresponding formulae in [KauOns] . 
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and Kaufman and Onsager observe that both Wn and Wn — t- as n — )■ oo, with Wn decaying at a 
much faster rate than Wn- They conclude that at T = Tc 

(100) (0-1,1 (Ti,i+n) ~ (-l)"cf W„ ^0 as n-^oo 

and so Mq = 0. In 1959, Fisher [Fislj picks up on these calculations, and using Stirling's formula 
in (j98p he shows via (jlOOp that as n — )• oo 

(101) a^,i+n) -Cn-^/\ T = n 

for some (undetermined) constant C. This is the first derivation of an explicit decay rate for the 
correlation function at Tc- Fisher argues further that off the lattice axes, (jlOip may be generalized 
(in the case Ji = J2) to 

(102) (ai,i a^+e,i+m) ~ 0(6) k-^/\ T = T, 

where tanO = £/m and k = (£^ + m^)^ — )• 00, and C{6) depends mildly on 6 (see more below). As 
Fisher learned later from Onsager himself (see reference 8 in [Fislj : also, cf. the reference preceding 
(j3ip above to [Dom2] ). the analogue of formula (jlOOp is in fact exact for correlations along the 
diagonal: 

(103) (ai,i cTi+„,i+„) = i-irWn, T = T,. 
Indeed from ([3T]) , ([32]) . for fcons = 1 at Tc, 

(104) (/.diag(e) = = e*^--')/' , < < 27r 
and so 

(105) ((^diag), = I , kez 

which agrees with (j97p (note from (j2ip that Wn is a Toeplitz determinant with shifted entries b^-i, 
and to observe (|103p one can write Wn as a determinant of the transposed matrix). If one repeats 
Fisher's calculation and keeps track of the constant, one obtains the exact asymptotics at T = Tc, 

(106) (c7i,i cri+„,i+„) = e3 A-^ 2T2 rT^I^ ~ ^ ^ + ' ' ') ^ n ^ 00 
where 

(107) A = eT2 e-^'(-i) 

is Glaisher's constant, (" = Riemann's zeta function. Formula ()106p was first obtained by T. T. Wu 
in |Wu) in serendipitous circumstances, as we describe below. The full expansion of ()106p to all 
orders was first given in [AuYPerl] . together with a similar result, in the symmetric case Ji = J2, 
for the next-to-diagonal correlation functions. 

For T > Tc, Fisher notes in |Fisl| that one can use the matrix viewpoint in [Onsl] to show that 
as k = \fW^rTn?' — )• cxd 

(108) (f7i,i ai+,,i+™) ~ B k-\ e-MT) 
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for some ^-dependent constant B = B{6), where to leading order as A; — oo 



(109) h{T)^kg{T). 

As T 4- Tc, g{T) \. 0. These calculations show, in particular, that Mq = for T > Tc- If one includes 
logarithmic correction terms in (|109p . one obtains (see |Fis2j ) for T > Tc 

(110) (^1,1 cjw+m) ~ c{e, T) (1 + 0(1)) 

as A: — )• oo. This is the first derivation of an explicit decay rate for the correlation function at T > Tc. 
In |KauUns] the authors show how the conclusion Mq = also follows from their determinantal 
formulae, at least for T > Tc sufficiently large. In |Kadj . Kadanoff obtains a decay rate for T > Tc 
of the same form as in (jllOp . but only in the double scaling limit T — T,, J, 0, A; — t- oo such that 
|r — Tel k = c= constant. In |Rya| , Ryazanov obtains the k~^^'^ decay rate at T = Tc, and he also 
obtains the general form (jllOp . at least in the case that T is large, T » Tc- The double scaling 
limit, T ^ Tc and k ^ oo, for the Ising model is an important issue, to which we will return at 
various points in this paper. The importance of this issue rests in the fact that many believe in 
universality for spin models in the sense that the double scaling limit for the spin-spin correlation 
function in the Ising model also gives precisely the same limit (modulo lattice-dependent factors) 
for a wide class of "non-integrable" two-dimensional models with short range interactions. 

Remark 4.2. The result that Mphys, the physically fundamental spontaneous magnetization defined 
above, is in fact zero for all T > Tc, was proven rigorously only many years later by Lebowitz in 
|Leb) in (1972). The proof in [Leb] uses a variety of techniques from statistical mechanics, but as 
shown in |LebM-Lof| |BGJ-LS| . Mphys = Mo = lim„_j.oo (f i,i 2 for all T > 0, and so the 

proof that Mphys = for T > Tc (in fact, for T > Tc) follows, a posteriori, from a proof that 
lim„^oo Dn{iPons) = 0, T > Tc. 

Remark 4.3. In the isotropic case, Ji = J2, explicit expressions for the correlation functions 
o"i+m,i+£) for values of {m,tj other than (0,^), (m,0), or {m,m), were derived by Shrock 
and Ghosh |ShGh] . In particular, they obtained expressions in the cases {m,i) = (2,1), (3,1), 
(3,2), (4, 1), (4,2), and (4,3), in terms of complete elliptic integrals K and E. They also inferred 
a general structural formula for arbitrary 1 cri_|_m,i+£) in terms of these elliptic integrals K and 
E. 

In further developments, Au-Yang, Jin, and Perk |AuYPer2] [AuYJFer] showed that the next- 
to-diagonal correlations (cJi^i cri^n,n) have a "bordered" Toeplitz determinant form, that is, an n x n 
Toeplitz determinant with the last column replaced by a certain explicit vector {bn-i,bn-2, ■ ■ ■ , bo)'^ ■ 
This representation is then used to express (c7i,i (7i+n,n) in terms of a solution to the Painleve VI 
equation by Witte in [Wij . 

5 New questions for SSLT 

In view of the above calculations for T > Tc, the Ising problem raised new questions and challenges 
in Toeplitz theory. From ([MI), (cji,i o-i,i+„) = -D„(v7ons) (similarly, from ([H]), (cri,i ai+„, = 
D„ (Vdiag))) ^-nd so we want to know precisely how L)„,((/?ons) behaves as n — t- 00 for T > Tc. But 
for T > Tc, we have < 71 < 1 < 72 (see (f27l) . (i30]l ). and so in contrast to the case T < Tc, the 
winding number of (/9ons is —1 7^ 0. And at T = Tc, we see from (f29|) that 9Jons(e*^) has a jump 
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discontinuity at ^ = 0. For such discontinuities the Fourier coefficients {logifoas)k decay as and 
so Yl'kLi ^ l(logv?ons)fcP = oo. In both cases SSLT, or more properly Theorem [71 does not apply. 
In 1966, in a tour de force in |Wuj . Wu gave the first derivation of the precise asymptotics of 
o"i,i+n) = Dni^ons) as n — )■ OO, for T > Tc- For T > Tc, Wu showed that for n — )■ oo 

{cTi,i ai,i+n) = (vrn)-5 (1 - 7?)^ (l " Ta"')"^ (1 " 7i72)"^ 

(111) A ^ ^1 ^2 ^ A3 ^ 

\ n n-^ n-^ 

with explicit expressions for Ai , A2 and in terms of 71 = zi z\ and 72 = -z^ /-^i > defined in (j26p . 
As T > Tc <;=^ 72 > 1, (jllip agrees in the order of decay with (jllOp for £ = and m = n, but now 
we have the precise rate and constants. For T = Tc, Wu obtained 

^^^2) (^I'l "^1.1+") = '^"^ (1 + 7i)^ (1 - 7i)~^ 

X (l + 5in"2 + 0(n~3)) 

with an explicit expression for Bi in terms of 71, and A is again Glaisher's constant ()107p . For 
T < Tc, Wu also obtained higher order terms in (pjj) . i.e., as n — )• 00 



(113) = (1 - /cL) M 1 + (27r n^) ' 7I" (72"' - 72 



Cl Co 

1 + - + 4 + 



with explicit expression for ci and C2 in terms of 71 and 72. In particular as 72 < 1, we see that 

o"i,i+n) ~^ (1 ~ ^ons) * — ^0 exponentially fast. 
There is an interesting twist to the story of Wu's computation of ()112p at Tc. As a warm-up 
problem at Tc, sinh j^J^ sinh = 1, he considered the limiting case when Ji — )• and J2 — )• oo 
in such a way that 71 i and 72 t 1 in (125P . and hence 



^i{n-e)/2 < 6) < 27r. 



1 _eie 
Wu recognized det ( 

ons ) j — k I ^ 

Cauchy determinant and computed the asymptotics 

V / 0<j k <n-l 

which we have recorded above in (jl06p . Wu was unaware at the time, just like Fisher before 
him, that he had serendipitously computed the correlation function along the diagonal at Tel Had 
Kaufman and Onsager published their formula for {ai^i ai+n,i+n) , this comedy of errors would 
probably have been avoided. 

In the symmetric situation, Ji = J2, one finds that at Tc, 71 = 3 — 2\/2 and so ^ir^^^ = 2*. 
Substituting this relation into (jll2p . we obtain as n — )• 00 



(114) (cTi,i C7i,i+„) = e3 2^/24 A-3n-3 (^1 + 0(^-^2) j, Ji = J2. 
This relation and (|106p are consistent with the formula (cf. (jl02p ) 

(115) = 6^2^/24^-3^-3 (^1 + (^^^^ , Ji = J2 
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where k = \J P ^ m? — )■ cxd. This asymptotic formula was derived by Cheng and Wu [ChenWu] in 
1967 (they also consider Ji / J2), but it was only established rigorously in the last year by H. Pinson 
[Pinj using a Phragmen-Lindelof type argument to interpolate between (|106p and (|114p . Recently 
two additional independent proofs of the rotational symmetry of (cii^i cri+^^i+m) as + — )• 

CO have been given by Dubedat [Dub] and by Chelkak, Hongler, and Izyurov |CheHonIzy| . In [Dubj 
the author proceeds by equating the product of two Ising correlators with a free field (bosonic) 
correlator. In |CheHonIzy| , the proof is based on convergence results for discrete holomorphic 
spinor variables pioneered by S. Smirnov. In earlier work, Smirnov used such variables to prove 
conformal invariance and to detect Schramm-Loewner evolution within the Ising model, in an 
appropriate scaling limit (for references, see |CheHonIzy| ). If Ji 7^ J2, one expects elliptical rather 
than rotational symmetry in the analog of (jllSp (see [ChenWu] and also (|199p below, which is 
taken from [WMTB] ). 

The full expansion of (cii^i ai+^^i+m) as £^ + m? cxd at T 7^ Tc in terms of so-called form 
factors f^\, is discussed in detail in [McCl] . 

(116) (fTi,i ai+m,i+e) = (1 - 0^ (^1 + f; f!n/^ for T < T, 
and 

00 

(117) (^Ti,i^i+m,i+^) = (l-t)3 5]/i'r'^ for T>T,. 

n=0 

For T < Tc, t = A;o„s and for T > Tc, t = k^^^. These expansions were first obtained by Wu, McCoy, 
Tracy and Barouch [WMTB] in 1976, but exphcit expressions for the form factors were only given 
in 1999 and 2000 by Nickel [Nickl] [Nick2j . 

The derivation of the form factor expansions (jll6p (|117p for (di^i ai+rn,i+i) was obtained in 
[WMTB] using the Fredholm determinant representations of Cheng and Wu [ChenWu] . In the 
special cases m = 0, ^— t-oo, orm = £— )-oo, Lyborg and McCoy [LyMcC] used the finite Toeplitz 
determinant representations to derive form factor expansions as in (jll6p (jll7p . but the detailed 
form of the expansions they obtained differs from the general expansions obtained in [WMTB] . It 
is an interesting, and still unresolved, question to reconcile the form factor expansions obtained in 
these two different ways. For recent results on the form factor expansion of the diagonal correlations 
cTi+m.i+m) which utilize results from [BorOk] and [GerCase] . see [WiFo] . 

6 Fisher-Hartwig singularities. Fisher-Hartwig conjecture 

In addition to the new questions and challenges raised in Toeplitz theory by the Ising model for 
T > Tc, other kinds of singularities also began to appear. For example, in 1964, in his work on the 
ground state of the one-dimensional system of impenetrable bosons [Lenl] Lenard expressed the 
one-particle density matrix for the system as a Toeplitz matrix with a symbol (/^(e*^) that vanished 
at some points on (see (1125P and particularly Remark 16 . 1 1 b elow for further discussion). In 1968, 
in a major step in the development and generalization of SSLT, Fisher and Hartwig [FisHartT] 
introduced a class of singular symbols for Toeplitz determinants, that allowed for zeros, (integrable) 
singularities, discontinuities and non-zero winding numbers. 
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The symbols of Fisher-Hartwig (FH) class have the following form (we use the notation in 
[DIKIQ 



(118) f{z) = e^(^) z^^-o ft- Y[\z- Zj\^^^ {z) z] 

j=0 



■ft 



z = e 



< 61 < 27r, 



for some m = 0, 1, 2, . . . , where 

2i = e'^^ J = 0,1, 



(119) 
(120) 

(121) 



= 9o<9l<--- <9m<27T, 

< arg z < 6j , 
9j < arg z < 27r, 



1 



^"i > "2' /^i ^ i = 0, l,...,m. 



and y(e*^) is a sufficiently smooth function on 5^. Here the condition on 5?aj insures integrability. 
Note that a Fisher-Hartwig singularity at zj, j = 1, . . . ,m, consists of a root-type singularity 



(122) 



I |2a, 
\Z-Zj\ ' 



2 sin ■ 



2a, 



and a jump singularity zft' gj^. (z) at (observe that zft g^j^. (z) is continuous at z = 1 for j 7^ 0). For 
j = 0, zq = 1, in addition to the root-type singularity (|122p at 9o = 0, one has 5^0, /3o(-^) ~ e~*'^'^°, 
< 61 < 27r and so z^" 5'^o,ft(^) = e* , < 6* < 27r, has a jump at z = zq = 1- The factors 
Zj are singled out to simplify comparison with the existing literature, in particular, [FisHartT] . 
Here are some examples of FH symbols: 

• If f{e^^) = -i for < 6* < TT and +i for tt < < 27r (cf [BasTr] ). then zq = 1, zi = -1 = e'"" , 
qq = «! = 0, /3o = ^ and /3i = — |, and 



(123) 



= 9i l{z) g^i^iiz) e-/2 



If /(e*^) is a sufficiently smooth function on (in particular, a trigonometric polynomial) 
with two zeros at < 6^1 < 6^2 < 27r (cf |BGrM] . |DIK2] ). then with zi = e^^\ Z2 = e''^^ , 
ai = 02 = 1/2, /3i = i, /32 = -i. 



(124) 



/(z) = e^(^) z''^+'3^ 



|z - zil |z - Z2I i[z) g 



iw(ii 

z - zi| |z - Z2I 1(2:) 9z2,-\i^) (^1/^2)^2 



for a suitable function y(e*^) on . 

The following symbol arose in Lenard's work in [Lenl] on impenetrable bosons mentioned 
above. 



(125) 



z — e 



i9l 



z — e 
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where 9i / 6*2 (mod 27r). Here zi = e*^i, Z2 = e"^^, ai = Q2 = 1/2 and /3i = /32 = 



26*2 
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• For the Ising model at T < T^, the winding number of fons (e*^) is zero and so (pons{z) = e^^^^ 
for a suitable smooth function y(e*^) on S^. For T = Tc, we have from (j29p with < 71 < 1, 

where = 1, ao = 0, /3o = — ^ and y(e*^) is smooth on S^. For T > Tc, we have from (f28]l 
with < 71 < 1 < 72, 

^ f.-^l (l-7ie-^)(72-e-^) / , 
(127) ^°"'V yl - |l_^,e-^| 11-72 e^'^l V ^ 

where = 1, ao = 0, /3o = —1 and again V{z) is smooth on S*-*^. 

From the point of view of Toeplitz determinants with FH singularities, Wu was the first to obtain 
|Wuj detailed asymptotic results for symbols with discontinuities as in p26|) and for symbols with 
non-zero winding as in (jl27p . Lenard was the first to obtain [Lenlj |Len2j detailed asymptotic 
results for symbols which have zeros on as in p25|) . For symbols (|118|) with 13 j = and aj real, 

> — ^, j = 0, . . . , m, Lenard made a general conjecture about the asymptotic form of 

Dnif), 

(128) !?„(/) ~S(e^, ao,...,am, 60, . . . ,0^) n^T=o-^'j e^^° 

as n — )■ 00, where Vq = 2^/0 V^(e'^) d9 as in (i37|l (i87|l . and i^'(/, ao; • • • > ctm, Sq, ■ ■ ■ , dm) is some 
(unspecified) constant. Moreover he was able to verify this conjecture for the case /(e*^) = 
\z - e*^i| \z - e*^2| in (fT25]l . 

Remark 6.1. As noted above, the symbol (I125P arose in Lenard's study of impenetrable bosons 
[Lenlj . More precisely, Lenard considered N free bosons in one dimension confined to a box of 
length L with periodic boundary conditions and lying in the (symmetric) ground state of the 
system. According to the criterion of Penrose and Onsager [PenOns] . Bose-Einstein condensation 
takes place if the largest eigenvalue A^^^^ of the density matrix pn,l{x, x') for the system is of order 
N as N L ^ 00. For the system at hand pn,l is translation invariant, pn,l{x, x') = pn,l{x — x'), 
so that the eigenvalues of pn^l are just the Fourier coefficients p^^\ of pn,l{-)- A simple calculation 

shows that A^^^^ = p^^\, the zeroth Fourier coefficient. Writing pn,l{£,) = xRn{'^'^C/ L), we have 
that 

^(max) ^ 

(129) ^ = 

(Note that in this case A^^^'' is actually independent of L.) Now it turns out that 

(130) RNit) = DN-i{ft) 



^ Note that although [Len2) was pubUshed in 1972, a prehminary unpubhshed version of the paper was already in 
existence in 1968 (see [FisHartl] ). 
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where ft is precisely a symbol of type (jl25p . 



(131) ftie 



Ad 



(Formula (jl30p was also obtained independently by Dyson, see |Lenlj .) At that point (April 10, 
1963, according to a letter Lenard has kindly made available to the authors) Lenard asked Szego 
for help in evaluating the Toeplitz determinant -Djv-i(/t)- On July 10, three months later, Szego 
wrote back to Lenard with the information that, although he had not been able to obtain precise 
asymptotics, he had obtained the bound 



(132) \RN{t)\ < 



Substituting (fTB^I) into we see that 

(max^* 

A 



eN 



sin(t/2) 



1/2 



(133) = O ^0, as N^oo. 



N 

In a letter of thanks to Szego on July 23, Lenard describes the situation as follows: "This has 
great physical interest: It may be expressed as saying that for the model considered there is no 
Bose-Einstein condensation in momentum space" o 

Following a further suggestion contained in the above correspondence with Szego in 1963, Lenard 
considered, in particular, symbols of the form 

(134) /(A,/^)(ei^) = |e*«-e^^/2|A|e'^-e"*^^Y, A,^>0. 

He started with the general formula for positive symbols /(e*^) > 

n— 1 

-2 



(135) DM) = ll^k 



k=0 



where Xfc > is the leading coefficient of the orthonormal polynomial pk{z) = Xk z ~\- ■ ■ ■ , k > 0, 
with respect to the weight of dfi = f{e'^) d0/{2TT) on as in 1^. For /(^-/^)(e*^) in (fTMl) . the 
polynomials Pkiz) can be expressed in terms of classical Jacobi polynomials on the interval [—1,1]. 
This then leads to explicit formulae for the Xfc's in terms of Euler's gamma functions. 

After Lenard published the bound (I132p for the impenetrable boson problem (a Plasma Physics 
Laboratory preprint of [Lenlj appeared in 1963), Dyson |Dy3] took up the problem of computing 
the asymptotics of (I129P precisely as — )• oo. In order to control the singularity in the symbol 
/t(e*^) = |e*^ — e**/^||e*^ — e~**/^|, Dyson resorted to his interpretation of L'Ar_i(/t) as a Coulomb 
gas at inverse temperature /3 = 2 with N — 1 unit charges free to move around on the unit circle, 
but now with two additional half-charges. The probability of finding the half-charges at an angular 
distance t apart is then proportional to P{t) = |2sin ||^/'^i?Ar(t). On physical grounds, Dyson then 
argued that for distances t much greater than the "Debye length" A = 2tt/N, the two half-charges 
are effectively screened off from each other, and so P{t) = const, for t >> A. In other words. 



*This confirms an earlier result of Schultz |Scli] proving the absence of condensation. Schultz used other methods 
which produced the (weaker) bound 0{N~'^^^ ). 
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RN{t) = k]y (sin I) for t » N ^, where kN is a constant. However, clearly = RNi^^) 
which can be evaluated explicitly as — )• oo since 



/t(e*^)|t=^ = /("'^)(e^^)U= 



(see discussion of (1134p above). Substituting the result into (I129P and neglecting the contribution 
from |t| < N~^, Dyson obtained, as N ^ oo, 

(max) / \ 2 

(136) ^ ~ CN'/\ C = (^) 2-^''A-^T Q 

where A is again Glaisher's constant (jl07p . The crux of the above calculation is clearly the explicit 
evaluation of i?Ar(7r) as — t- oo. 

More detailed analysis of ()129p requires further study of the double-scaling limit as t — >• 1, 
N — )• oo. We return to this issue in Section [TTl 

Inspired by the work of Wu in |Wu] . and perhaps also influenced by the calculations of Lenard 
in [Lenl] |Len2] . Fisher and Hartwig [FisHartT] made a conjecture about the asymptotic form of 
Dn{f) for general symbols (IllSp . 

DnU) = E (e^ ^0, • • • , /3m, Gq,. . . , 9 m) 

^^^^^ ^oM-m e-Vo (1 + o{l)) , Vo = ^ fj V(e^') d9 



as n — )• cxD, together with a conjecture on the general form of the constant E. In |FisHart2] . Hartwig 
and Fisher present various examples, heuristics and theorems in support of their conjecture. 

In the remarkable paper |Wid2j . representing "a jump of several quanta in depth and sophis- 
tication" (see MathSciNet MR0331107), Widom verified Lenard's conjecture for complex aj with 
Kaj > — |, j = 0, . . . ,m, and y(e*^) suitably smooth, and obtained an explicit expression for the 
constant 



(138) 



E{e^,ao,...,am,eo,...,9m)=E{e^) J] - e''^ 

0<i<fe<m 



i=0 jr=0 



where 

(139) E{e^) = e^k=ikVkV.,, ^ ^ Fourier coefficient of v(e'^ 



(140) V\^e"'j=V\^e"'j-VQ 
and 

(141) Eo, = G'^il + a)/G(l + 2a), G{z) = Barnes G-function (see [Bar]). 

Formula (jl38p is consistent with and confirms the general form of the conjecture |FisHartT] for the 
constant term in the case /3o = /3i = • • • = f3m = 0. 
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In [Baslj Basor considered complex aj with ^aj > — ^, (3j pure imaginary, j = 0, . . . ,m, and 
y(e*^) suitably smooth, and verified (|137p with constant term, again consistent with the general 
conjecture in [FisHartl] . 



E [e ,ao, . . . ,am., f3o, . . . , (3m,0o 



J ■ ■ ■ 1 



Eie' 



j=0 



(142) 



X n 

0<j<k<m 



Zj Zfc 



Zi e' 



ajl3k-akl3j 



i=o 



G{l + aj + /3j)G{l + aj-pj) 



G{l + 2aj) 



with -E(e^) and G{z) as above, and 

(143) b+{z)=e^^=^^'^'\ 



E — oo 



The branches in (jl42p are determined naturally, b^{zj) 



Zi e' 



(ajl3k-akl3j) 



eC-^.+ftOE.-liV'.^ etc, and 



< 0j < 6*^ < 27r. 



In the paper |BasHelt] discussed above, Basor and Helton also introduce the first example of a 
so-called "separation" theorem of the following form: For symbols and ip of a certain prescribed 
type, including FH symbols with disjoint singularities, they show that 



(144) 



hm ^4^^ 



exists. Moreover, they provide an explicit form for L{(p,%l)). Iterating this result, they are able 
to reduce the problem for general FH symbols to adding in pure FH singularities of the form 
z^\z—z\^'^g'l^p{z)z~^, \z\ = 1, one at a time. They analyze the asymptotics of Toeplitz determinants 
with such pure (i.e., with 1/ = 0) FH singularities in the following way. In addition to the results 
described above, in |Wid2| Widom also verified (I137j) for a single FH singularity with |5}?q| < ^ and 
< i, but without determining the exact value of the constant E. But then Basor 's formula 



()142p for E obtained in [Baslj for 'Raj > — 2 and j3j pure imaginary, together with Vitali's theorem, 
directly yields the desired explicit asymptotics for a single pure FH singularity with |3f?a| < ^ and 
< |. (In fact, as discovered later by Bottcher and Silbermann, see ()145p below, there exists an 
exact formula for a Toeplitz determinant with a single pure FH singularity in terms of the Barnes 
G- functions, which easily yields the asymptotics). In the end, Basor and Helton are able to verify 
(|137p (I142P for Oj, 13k in the open set {maxj,fc(|aj|, |/3fc|) < 6} C C^" for some (small) 5 > 0. 

In [Bas3] . Basor showed using a separation theorem that (jl37p (|142p also hold in the case Oj = 0, 
|5R/3j| < i,i = 0,...,m. 

Bottcher and Silbermann [BottSilbS] finally proved (jl37p with Basor's constant (jl42p under the 
assumption that |Kaj| < 1/2 and < 1/2 for all j. This is the most general setting in which 

the conjecture ()137p is true in its original form, that is, with the exponent of n conjectured to be 
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^(oj — Pj), without the need to take care of degenerate situations and of the ambiguity in the /3's 
(see (|154p below). We note that in [BottSilbB] it had aheady been shown that (jl37p does not hold 
in general if aj+/3j and aj — (3j are nonnegative integers for all j = 0, . . . , m, and a corrected version 
of the conjecture was proved in this case (this is a particular case of a formula conjectured by Basor 
and Tracy: see next section). The developments up to the 1990's are described in [Bottl] and [Ehr] 
and also in the books [BottSilb2| |BottSilb3] . We note that certain separation theorems also play 
a key role in the work of Bottcher and Silbermann. Another important element in their proofs 
is an explicit formula, noted above, for a Toeplitz determinant with a single, pure FH singularity 
found by the authors in [BottSilbS] (see also [BottWid2] [BasChej for later alternative derivations) . 
Namely, if V{z) = 0, m = 0, then with oq = a, /3o = /3, 

G(l + Q + /3)G(1 + a - /3) G(n + l)G(n + l + 2a) 

^n\j) — 



^-^45) G(l + 2a) G(n + l + a + /3)G(n + l + a-/3)' 

3fJa>-^, Q±/3 / -1,-2,..., n>l. 

The condition on a it /3 is needed because the G- function has zeros at z = 0,-1,-2,... The 
asymptotics for (|145p follow from the asymptotics of the G- function (see [Bar] ) 

logG(t + a + 1) = - log^ - ^ - at + log(27r) 

(1^6) ft- a' l\ 

+ ( ^ + at + y - — j logt + o(l), ast-^oo 

where again A is Glaisher's constant ()107p . The resulting asymptotics are consistent with (jl37p 
(fn2]) . 

In view of the Heine representation (j45p . it is not surprising that the formula (|145p is related 
to the Selberg integral. In fact, independently of BottSilb^ . a more general formula was obtained 



by Selberg himself: 

f27r /■27r 

27r 



f 

Jo 



'0 Jo 

(147) 



/■^i" . . i2'v - : 

I n h-^'l n 

r(l + 2a + j7)r(l + (i + l)7) 



2" d9j 



nfa 



j=o - - + " + + ^^^^^^ + a - /3 + n)r{i + 7) 



where 



^ 1 ^ . fl K(2a + 1) 

> --, K7>-mm<^-, 



2 I n n — 1 



In a published form, this formula first appeared in the thesis of Morris in 1982 (see [ForWar] 
for a historic account). The interest in ()147p was due to the efforts to prove one of Dyson's 



conjectures in the theory of random matrices in Dyl . If we set 7 = 1 in (11470 . and use the relation 
G{z + 1) = T{z) G{z), we recover (fHS]) from (liSll . 

The Barnes' function G{z) is a basic object in analysis, bearing the same relation to the gamma 
function as the gamma function bears to the function z, T{z + 1) = zT{z). Just as logr(z) is a 
sum of values of logz, so logG(2;) is the iterated sum. It is of interest to delineate some of the 
analytic pathways by which G{z) appears in the asymptotics of Dn{f). 
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The most direct example, a particular case of ()145p . arises for the diagonal correlation function 
at T,, where (^diag(^) = e^("-^)/2, i.e. V = 0, m = 0, ao = 0, = -1/2, in 
gives rise to the Cauchy determinant (jl03p which can be evaluated as in (j98]) in terms of a product 
of gamma functions, and hence in terms of G-functions, ng=i ^('7 + 1) = ^(n + 1)/G(2), etc. 

For a system of polynomials orthonormal with respect to a function / on the unit circle, the 
leading coefhcients Xk are related to Dn{f) by the product formula (llSSp . As we discussed above, 
in the case of / given by (jl34p the coefficients Xk are expressed in terms of gamma functions, 
and therefore G-functions appear in the product. More generally, there is a certain universality in 
the asymptotic behavior of orthogonal polynomials associated with FH symbols on (see |DIK1] 
[DIKSj ). and, mutatis mutandis, for general FH symbols /, the G-functions enter the arena through 
the same stage door. 

As noted above in (1126p . at T = Tc, (/?ons(e*^) has a single FH singularity at zq = 1 with 
ao = and /3o = —\- The first confirmation of (jl37p (jl42p which allowed for these values was 
given in |BottSilb4] . where the authors proved the result for m = 0, 5Rao > 0, 3ft (ao + /3o) > 
— 1 and 3ft(ao — /3o) > —1- It is instructive to recover Wu's asymptotic formula (jll2p from 
(fT37D (fTi2]l . Here V{e'^) = \\og[j^^^^ and we find = and Vk = -if^ /2k for k / 

0. Thus E{e^) = 6-3^"=! '=^1"''/'=' = (l-7i2)2. On the other hand (6+(l))"5 (5_(i))5 = 

e3Er=i7fAeiSfc"-i^i"/l'=l = (l-^^)-!. But G {\) = 221 7r-V4 ^-3/2^ ^j^^^g ^ ig 

Glaisher's constant (see jBarj ). and we have G (|) = T (|) G (i) = 2m vr^ Thus 

G (i) G (I) = 21^ e3 A-'^. Substituting these relations into (fTl2D . we obtain £;(e^, 0, i, 6^0 = O) = 

2 12 64 A" (1 + 71)4 (1 — 71) 4, and as n~"'-'> = n~i, we arrive at precisely the leading order term 
in (fTT2D . 

The definitive result for the leading order asymptotics of Toeplitz determinants with FH 
singularities was obtained by Ehrhardt in his PhD thesis at TU Chemnitz in 1997 (see [Ehr] ) . Let 
denote the seminorm 



(148) III/3III =max|3ft/3,--3f?/3fc| 

j,k 

where 1 < j. A; < m if oq = /So = 0, and < j, k < m otherwise. If m = 0, set |||/3||| = 0. 
Note that in the case of a single sing ularity, we always have |||/3||| = 0. If V{e^^) is G°° on S\ 
|||/3||| < 1, 3ftafc > — ^ and aj it f3j / —1, —2, ... for j. A; = 0, 1, . . . , m, then Ehrhardt proved the 
Fisher-Hartwig conjecture ()137p with Basor's constant (|142p . As above, the conditions on aj it f3j 
ensure that Basor's constant is non-zero and (jl42p indeed provides the leading order asymptotics 
for Dn{f) as n — )• 00. A key element in Ehrhardt's proof is a suitably generalized version of the 
separation theorem in |BasHelt| FI 

An independent proof of Ehrhardt's result was given recently in |DIK3j where only a finite degree 
of smoothness is needed for y(e*^). In [DIK3| . estimates for the error term in the asymptotics are 
also provided. If V{z) £ C°° on S*^, the asymptotics (|137p (jl42p hold with error term 



(149) 0(1) = 0{n 

Note that the smallness of the error term breaks down if |||/3||| > 1. We will consider this situation 
in detail in the next section. The main steps of the method of jDIK3j are as follows. First, the 



^Note that the RHS of (fT57)) and (fTi^ make sense for any a^, 2qj / -1, -2, . . . . In fact, ((157)1 (fTi^ still hold in 
this case. This extension to the left of the lines Kqj = —1/2 is explained in Ehr and corresponds to replacing the 
symbol / by an appropriate distribution. 
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authors derive differential identities for tlie logarithmic derivatives of Dn{f), in particular, for 

(150) A log !)„(/), ^logDM), 

oaj opj 

in terms of the associated OPUC's of degrees n and n + 1 only. Second, the authors obtain the 
asymptotics of these OPUC's using their Riemann-Hilbert representation and the steepest descent 
analysis of the Riemann-Hilbert problem for large n. Third, the asymptotic formulas so obtained 
are substituted into the RHS of the differential identities, and the final result ()137p ()142|) is obtained 
by integration In |DIK1| . [DIK3] the authors also note various uniformity properties of the large 
n asymptotics for Toeplitz determinants in the parameters, such as aj, Pj, on which / depends. 
This issue becomes important in discussions of double-scaling limits (cf Section [TT|) . 

For the Ising model at T > Tc, we see from ()127p that zq = 1, oq = and /3o = —1, and so 
ao + (3q = —1. Thus (|137p (|142p do not provide the leading order asymptotics for (cji^i (Ti_i-|_„) in 
this case. One must analyze lower order terms, in fact go "beyond all orders" in the language of 
Kruskal, to capture the exponential decay in Wu's result (jllip (see Remark l7.ll below). 



7 Basor- Tracy conjecture 

As already observed in [BottSilbS] |BottSilb6] . Ehrhardt's condition on the seminorm, |||/3||| < 1, is 
not just a technical matter. In a critical calculation in 1991, Basor and Tracy [BasTrj considered a 
particular symbol /(e*^) which has two FH singularities at = 1 and zi = e*'^ = — 1 respectively 



(151) f[e'')=9,,l{z)g_,^^^{z)e^-/^ 

with /3o = ^, /3i = — ^ so that |||/3||| = 1. By direct computation they found that as n — oo 



1 + (-1)" 12 ^ flV ^ /3^ ^ 



(152) DM) = ^^^-G^-) Gl^-j (1+0(1)) 

which is not of the general FH asymptotic form. Basor and Tracy observed, however, that /(e'^) 
had another FH representation with /3o = — ^ and /3i = i, 

'in/2 



(153) f[e'")=e'^g,_^{z)g_,^{z)e 

If one substitutes the values F = 0, oq = «i = 0, /3o = ^, /3i = — ^ into ()137p (jl42p . one obtains 
for dm]) 

n-ii+i) |i-e-| ^ G(l) G(l) + 

= (2n)-^G(|)' G(i)'(l + o(l)) 

as G(l) = 1. On the other hand, if one substitutes the values V = iT:, = ai = 0, = — i, 
/?! = |, one obtains for ()153p 

e--(2n)-^G2(l)G2(|)(l + o(l)). 



^'^ Note that if one just used the formula p35|) instead of the differential identities, then the product Hfelo Xfe ^ 
small no, and hence the constant factor E in (|137[) . would remain undetermined. 
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Just adding these two asymptotic forms, we obtain precisely the Basor- Tracy result (jl52p . So we 
see that if a symbol /(e*^) has more than one FH representation, then the asymptotics of Dn{f ) 
is given, at least in this case, by some combination of the Fisher-Hartwig asymptotics for each of 
them. Based on this example, Basor and Tracy [BasTrj made the following bold, general conjecture. 
First note the following. Let /(e*^) be a Fisher-Hartwig symbol as in (jllSp . Let nQ,ni, . . . ,77,^ be 

a set of integers with YlY=o ""i ~ ^' f ('^' denote the FH symbol obtained by replacing /3j 

with = (3j+nj, 0<j<m, and e^(^) by (n"=o-2f ) = e^+'^T=o ".e.. Then f0) = f {z-p 
gives another FH representation for f{z)^ 



j=0 



. 



(154) f{z) = fi^z;(3), 

Moreover, it is easy to see that all FH representations of f{z) arise in this way. Given /3, we call 



3=0 



(155) Op -- 
the orbit of /3. We consider the discrete minimization problem 

(156) = min 

It is clear that the minimum is indeed obtained in O^. Let M 








0, 



. . . , m and all /3 G 7W ^ . 



We say A^/j is non- degenerate if it /3j 7^ — 1, — 2, . . . for all j 

Conjecture 8. [BasTrj Let f{e'^) be given as in (jllSp with y(e*^) sufficiently smooth on and 
> < j < m. Suppose is non-degenerate. Then as n — )• 00, 



(157) 



Dn{f)= J2 [i?n (/(^)) (1 + 0(1)) 



Each Rn yfif3)j stands for the FH asymptotic form (|137p (|142p with /3 and (^njLo ^^j place 
of /3 and e^ respectively. 

Note that by the definition of A4/3, all the terms Rn (/(/?)) in (|157p have the same order of 

and no+n-i = 0, 

2 



magnitude. Indeed for all ^ E M/^, EjLo =Pp- TJjU (^/^j) • 

To see how this conjecture works for /(e*^) in (jl5ip . where = \ and j3i 
consider 

1 , \2/-i \2 

j=0 

which clearly achieves its minimum at no = 



+ 



-\-no 



2 ( ^ + ^0 



correspond to /3 = /3 = (i, — i) and /3 = (^ — 1, — ^ + l) = (— ^, ^), respectively. Thus the leading 



asymptotics of Dn{f) is a sum of contributions from f{z 
(fT52D . 

Conjecture [8] was proved recently in [DIKl] . 



ni and also at no = — 1 = — ni. These values 
/(z;/?) and /(z) = / ( z; /3 ) , as in 
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Theorem 9. IDIKlf Conjecture\E holds. 

The relation of the semi-norm |||/3||| to the minimization problem is given by the following 
elementary result. 

Lemma 10. jPIKl^ Given ji, there are only two, mutually exclusive, possibilities: Either 

(i) There exists /3 G such that |||/3||| < 1- Then such a (3 is unique and it is the unique 
element of = |/3| 

or 

(a) There exists (3 G Op such that |||/3||| = 1. Then there are at least two such (3 's and all of them 
are obtained from each other by a repeated application of the following rule: add 1 to a f3j 
with the smallest real part, ^f3j = min^ K/3fc, and subtract 1 from a 13 ji with the largest real 
part, ?fi(3j> = maxfc 0k- Moreover = |/3 G O/3 : |||;5||| = l|. 

Clearly for /(e^^) in (USQ), |p||| = 1, if and only if;5 = /3or/3 = /3 + (-1, 1). Also we see from 
(i) that the case |||/3||| < 1 in Ehrhardt's Theorem corresponds to a unique minimizer in ()156p . and 
hence to a single term in ()157p . 

The proof of Conjecture [8] in [DIKlj proceeds as follows. Given /3, if there exists f3 G Op with 

|||/3||| < 1, then /3 is unique and as n — )• 00, Dn{f) = Rn {j{f^)^ (1 + (Ehrhardt's Theorem). 

On the other hand, if there exists /3 G with |||^||| = 1, then Mp = |^ G O/3 : |||^||| = l| has at 
least two elements. Then for some 6, 

(158) h<0j<b+l 

for all (3j, and for some p, ^ > there are p values of j such that = b and i values for which 
3^/3j = 6 + 1. Let /(e*^) be the FH symbol {not a FH-representation of /(e*^)) with Oj, V as before 
but (3j = /3j if ?fi/3j < 6 + 1 and = /3j - 1 if 0j = b + l. Then clearly 

(159) fiz) = c/fiz) 
for some simple constant c, and by a general relation 

(160, z,„(/,.£!tii;^ 

where 

(161) F„ = det (0)) 

V / 0<j,k<l-l 

and TTq^z) = z'^ + . . . ,q > 0, are the monic orthogonal polynomials with respect to the (in general 
non-real) weight /(e*^) ^ on 5^, 

(162) r (e^^) e-^ / (e^^) ^ = 0, < r < q. 
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(A part of the proof of Conjecture [8] consists in showing that the tt^'s exist and are unique for q 
sufficiently large, q > qj.) Now clearly |||/3||| < 1, and so we can evaluate Dn{f) as n — >■ oo using 
Ehrhardt's Theorem. The proof of the conjecture then reduces by ()160p (I16ip to evaluating the 
polynomials 7rn{z), . . . , TTn+e-i{z) and their derivatives at z = 0, as n — )• oo. This is effected by using 
the steepest-descent method of Deift and Zhou for Riemann-Hilbert problems mentioned in Section 
[3] above (see also, for example, |DKMVZT] . |DKMVZ 2^). The fact that orthogonal polynomials 
with respect to a weight on the line, can be expressed in terms of a Riemann-Hilbert problem is 
due to Fokas, Its and Kitaev in [FIK] : the extension of this result to orthogonal polynomials on 
the circle is given in |BD J] . 

Remark 7.1. As noted above, for the Ising model, at T > T^, 93ons(e*^) in (|127p corresponds to 
a degenerate case for Ehrhardt's Theorem. However (/7ons(e*^) = — e~*^ (^(e*^) where ^(e*^) = 

'^^ '^ie il 1"^^ ^ ie| smooth, nonzero, and has no winding on S^. Hence Dn{(p) can be computed as 
n — )• oo using SSLT as in ([85]) . Analogous to (jl59p . (jl60p . we have the formula 

(163) Dn{iPons) = vfn(O) -Dn(^) 

where TTniz) = z" -|- ...,n > are the complementary orthogonal polynomials to the 7r„'s (see 
[DIKlj ) 

(164) £^ 5f„ (e-^) e^^' ^(e^^) ^ = 0, < j < n. 

Applying the steepest descent method to the Riemann-Hilbert Problem for vf^, and taking into 
account the analyticity of if{z) in an annulus around {|z| = 1}, we obtain the precise decay 
rate of vf„(0) as n — t- cxd, and hence via (jl63p . the precise leading order exponential decay rate 
gj^-i/2^-n _|_ Q^^i'j'^ [j^ (jllip . The steepest descent method also yields the higher order terms in 
the expression, to any desired order. An alternative, and very elegant, proof of (jllip using a contour 
deformation argument, is given in [FF] . For (/3diag(e*^) and the diagonal correlation function, the 
situation is similar, and one obtains: 



h-n rz 

(165) cTi+.,i+„) = ^ 1.-2, + T > T,. 

We also note that in Section 10.7 of jBottSilb3] and in the papers [BottWidT] |Wid7| . the authors 
also analyze the Toeplitz determinants arising in the case T > Tc- 

Remark 7.2. Formulae ()160p and ()163p should be compared respectively with formulae pOSp and 
(1971) in [FisHartlj . 



8 Hankel and Toeplitz+Hankel determinants 

In [DIKlj the authors also consider Hankel and Toeplitz + Hankel determinants with FH-type 
singularities. For a weight w = w{x) > on [—1, 1], the associated Hankel determinant is given by 
[5i5] 



(166) D^^) (w) = det ( / x^+'' w{x) dx ) 
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The Hankel determinant Dn (w) is related to the polynomials orthogonal with respect to the 
weight w on the interval [—1, 1] in a similar way as a Toeplitz determinant is related to OPUC's. 
For fixed r = 0, 1, 2, . . . we consider w{x) of the following form: 

r+l 

(167) w{x) = e^(^) n 1^ ~ '^i(^) 

3=0 

where 



1 — Aq > Ai > • • • > Aj.-|_i — — 1, 

—1 < X < Xj, 

l>X>\j 



LOj{x) 



= = 0, 



'Raj > --, 
2 



j = 0,l,...,r + l 



and where U{x) is a sufficiently smooth function on [—1, 1]. 

Remark 8.1. Note that there is no loss of generality in setting /3o = /3-r+i = and G 
as the functions uJo{x),ujr+i{x) are just constant on (—1, 1) and ujj {x; /3j + kj) = (—1)^^ ujj {x;(5j) 
for kj G Z. 

In |DIK1] the authors prove that for w as in (fT67|) with G (-i, i) , j = 1, . . . , r, 
(168) D[,^\w) = 2-"' G^n-i+2("g+-?)+E-=i("?-/3|) (i + ^(i)) 



,1 11 

2' 2\ 



as n — )• oo, 



for certain explicit constants Gh and i?/;- depending on U{x), Xj and /3j (see (1.37) in [DIK1| ). 
This result is derived by introducing the function 

(169) /(^e*^) = /(^e^^^'^^^)) ='w(cos6')|sine|, < 6 < Itt 

on S^. With w as in ()167p . / has FH singularities of the form 



cen<^")|z-l|4°o+i|^ + l|4a,.+,+i -Q 



\Z — Zj\ ^ \Z — Zj^ 



(170) 



where 



ll[9^,,.,,{z)z^^ 9,^^,^iz){z'^y'^ 



< 6* < 27r 



cos 



Aj, = 



^ i = 0,1,..., r + l, = ^0 < ^1 < ••• < 0r+i =7r 
j = 0,l,...,r + l 



and 
(171) 



Vie'") = U{cose). 
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The Hankel determinant Z)„ \w) and the Toephtz determinant D2n{f) are related in the fohowing 
way: 

4(."- vr2n(l) vr2„(-lj 

where 7rfc(2) = + . . . ,k > 0, are the monic orthogonal polynomials with respect to the weight 
/(e*«) ^ on as in (fT62]) . 

As the semi-norm |||/3||| < 1 for fie^^^, we may use Ehrhardt's result to evaluate D2n{f), n — )■ oo, 
and the result (jl68p follows from (jl72p by computing the asymptotics for 'K2n{z) for z = 0,ibl. 
Note that if 3f?/3j = \ for some j G {0, 1, . . . , r + 1}, then |3f?/3j - 3f? = 1 and so |||/3||| = 1. In 

this case we would have had to use the Basor- Tracy form p57|) for D2n{f)'- we have restricted our 
attention to the case G o^^Y foi" simplicity. 

Remark 8.2. For a discussion of results of various authors, particularly Basor and Ehrhardt, related 
to (fT68]l . see [raKlJ, Remark 1.24. 

For even FH symbols /, i.e., /(e*^) = /(e*^^'^^^)) , the authors in |DIKlj also consider Toeplitz 
+ Hankel determinants of four types defined in terms of the Fourier coefficients of /, 

(173) det [fj-k + fj+kT~k=o ' "^^^ ^fj-k ~ h+k+2Tj~k=Q ) det ± fj+k+iTj'k=o ■ 

Such matrices arise, in particular, in the theory of classical groups and its application to random 
matrix theory and statistical mechanics, and there are simple relations between these determinants 
and Hankel determinants on [—1, 1] with added singularities at the end-points (see [DIKlj and the 
references therein). For example, for /(e*^) = / (e*^^'^"^)), 

nn2-2n+2 

(174) det (/,„fc + /,+.);-,Lo = ^ D^n^ (^) 

where d'^\v) is the Hankel determinant with symbol v{x) = /(e*^^^^) /\/l — x^, 0{x) = arccosx, 
on [—1, 1], with similar formulae for the other three determinants. Let /[e"^^^ = f (e*^^'^"^)) with 
FH singularities at = < < ■ ■ ■ < < ^r+i = ^ and complementary singularities at 27r — 6j, 
1 < J < "T. Suppose that 3?/3j G {—^: ^); j = 1; • • • > and /3o = Pr+i = 0. Then utilizing ()174p and 
()168p . the authors show that as n — t- oo 

(^75) det (/,_, + f,+k)lklo = G^^^n^Ui-^-f^lhU-l+^h.'-o-ar^^) 

X Et+h {1 + o{l)) 

for certain explicit constants Gt+h, Et+h, with similar results for the other three Toeplitz+Hankel 
determinants in ()173p . 

Remark 8.3. For a discussion of the related results of Basor and Ehrhardt for determinants of type 
(fT73]) . see [DIKlj . Remark 1.27. 

Remark 8.4. In [FFj the authors discuss a large number of conjectures and results for Toeplitz, 
Hankel, and Toeplitz+Hankel determinants with FH-type singularities. The determinants arise in 
turn from problems in statistical physics. Many of the conjectures in |FF] can now be resolved 
using results from [DIKlj such as (fTHHl) and (fTTHl) . 
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9 Some applications 



As advertised at the outset, the goal of this paper has been to show how SSLT developed and was 
generalized in response to questions arising in the analysis of the Ising model. Although the case 
III/3III = 1 in [DIKl] does not arise in the Ising model (the same is true for the Hankel and Toeplitz 
+ Hankel determinants (jl66p and ()173p respectively) , this case does arise in many other problems in 
statistical physics. For example, the probability PEin) of a string of length n of ferromagnetically 
aligned spins in the antiferromagnetic ground state in the XY spin chain, is given, for a certain 
range of parameters, by a Toeplitz determinant with 2 /3-singularities such that |||/3||| = 1. Thus 
(jl57p verifies the result on PE{n), n — )• (X), presented in [Fr Ab| . which the authors based on the 
Basor- Tracy conjecture. In a similar way, the results in |DIK1] can be used to justify the asymptotic 
results for correlations arising in the theory of the impenetrable Bose gas that were obtained in 
the work of Ovchinnikov [Ov] on the basis of the Basor- Tracy conjecture. In another direction in 
[GGMj ■ the authors use ()157p to evaluate the long-time behavior of a variety of non-equilibrium 
ID many-body problems. In particular, in the Fermi edge singularity and tunneling spectroscopy 
problems considered in |GGM| . the authors show that the various contributions to (|157p have a 
very transparent physical meaning: they correspond to contributions to the Green's functions for 
these problems from multiple Fermi edges. 

Theorem [9] can also be used to analyze the asymptotic behavior of the eigenvalues a["^ < Ag"^ < 
• • • < An"\ n — )• oo, of a Toeplitz matrix Tn{'f) with a real symbol 99. At the beginning of our story 
(see Theorem [2]) we interpreted Szego's early results as saying that the A^""^ 's were equidistributed 

(n) 

as n — )• 00. But what about the behavior of individual eigenvalues, A^ , n — )■ cxd? Early work on 
this question considered the so-called extreme eigenvalues a[."^ and A^"^^ with k fixed as n — ?■ 00 

(see [GreSzj . [BGrMj and the references therein). The important case for A^"^ with n — )• 00 and 
A:/n — )• a: G (0, 1) was only considered for the first time recently in [BGrM| . Recall that by general 
principles (see e.g. |GreSz| ). for bounded, real- valued functions /(e*^), the eigenvalues {A^*^^} lie 
in the interval [L,M], L = inf /(e*^), M = sup /(e*"), and as n — ;> 00 the A^"''''s fih out the 
interval in the sense that if A G [L,M], then there exist kn{\) such that ~^ ^- I^i [BGrMj 

the authors consider function /(e*^) that are unimodal, real- valued trigonometric polynomials, i.e., 
for some < 6*0 < 27r, ^ f{e'^) > for < 6* < 6*0 and ^ /(e'^) < for 0o < f < 27r, with the 
additional property that ^ /(e*^) 7^ at = and 9 = 6q. For any A G (L, M), there exist unique 
< 6'i(A) < 6*0 < 6'2(A) < 2tt such that /(e^'^i(^)) - A = 0, j = 1, 2. The function 

extends to a homeomorphism from [L^M] onto [0,7r]. In [BGrMj the authors prove, in particular, 
that for any x G (0, 1) 

(176) Ai") = A. + 0(1) 

as n — )• CXD and k/n ^ where ip{\x) = ttx, and 

("^) = 

if \k/n — x\ = O (^) as n — )• 00. 
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The estimate ()177p brings precision to the equidistribution result in Theorem [2j In [DIK2] the 
authors extend the result in |BGrMj to C°° unimodal, real- valued functions on S^, and they also 
consider piecewise-constant symbols of the form 



(178) 



-id 



0,^1 u 



= 2777 



where Q < 6i < 62 < Itt and 7 > are given. For symbols such as (jl78p . the eigenvalues A^"^ 
cannot be spaced 1/n apart as in ()177p . This is clear, for example, from ([8]) in Theorem [21 if we let 
-F(A) be any non-negative continuous function with support in (l,e^'^'''). And indeed, the general 
result in |Bas2| asserts that in the case (jl78p . for any small e > 0, the interval (l + e,e^'"'^ — e) 

contains order logn eigenvalues A^"'^ spaced at distance 1/logn apart, as n — )• 00: The bulk of the 
eigenvalues accumulate at 1 and at e^'^"^ . On the basis of numerical computations, the authors in 
[LeSh] (see also |LeSoSo| ) conjectured that in the case 

(179) 02-91= 27rp/q, p, q e Z, 0<p<q 

there is a "near periodicity" in the spectrum of Tn{f). In |DIK2] the authors prove this conjecture 

1 + e, e^-T-e), 



in the following sense: for n sufficiently large, if A^"^ in an eigenvalue of T„(/) in 



then there is an eigenvalue A 



{n+q) 



of Tri+q such that 



(n) 



A 



{n+q) 



< c{n log n) . Of course 



l/(n logn) << 1/logn, the spacing of the A^"'''s, k = 1, . . . ,n. 

The proof of the above results in |DIK2| rests on the simple fact that det {Tn{f) — A) = 
det (r„(/ — A)), so that the eigenvalue equation for Tn{f) is equivalent to the vanishing of the 
Toeplitz determinant Dn{f — X). Now if /(e*^) is a smooth, unimodal function on , then /(e*^) — A 



is a FH symbol with two singularities at zi 
(see (HMD) 



Z2 



ai = a2 



1 

2' 



/3i = -h 



(180) 



zi\ \z- Z2\g,^i{z) 9,^^_i{z) 



where y(e*^) is smooth on 5^. This is a case of a FH symbol with |||/3||| = 1. Therefore we need 
to use Theorem [9] to obtain the asymptotics for Dn{f — A). The set Mg here contains 2 elements, 
and the condition for an eigenvalue is the vanishing of the sum in (I157p : 

(181) DM - A) = Rn{f{f3i,/32) - A)(l + o(l)) + i?„(/(/3i -1,^2 + 1)- A)(l + o(l)) = 

The asymptotic results on the eigenvalues follow from this formula. 

The Toeplitz matrix T„(/) with the symbol /(e*^) in (|178p is handled similarly. The function 
f(z) — X in this case has 2 /3-type singularities whose positions on the circle are fixed, but the 
imaginary parts of /3's depend on A: /3i = — /32 = ^ + ^7(A). Note that we again have = 1, 

and the condition for the eigenvalues is thus given by Theorem [UJ 

Finally we describe an example where the results in |DIK1] are needed for a Toeplitz-I-Hankel 
determinant of type (|173p . The following problem arises in the framework of the random matrix 
approach to the theory of the Riemann zeta-function and other L-functions (see [Keatj for a recent 
survey). Define 



2k 



(182) 



J9 



2 sin ■ 



.V{e^O) 



A; e N 
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where 



Vie 



id 



2k 



u{y) < 



f2 Ci(|0 + 27rj|(logy)(logX)) Uy-log 



Ci(z) 



J = -0O 

°° cost 



2 sin ■ 



dt 



and u{y) is a smooth, non-negative function supported on 
Consider the fohowing average over the orthogonal group SO{2n), 



and of total mass one. 



(183) 



E 



SO{2n) 




Here e^*^^, . . . e^*^" are the eigenvalues of a random matrix in SO{2n): note that 4> is even, </>(e*^) = 
0(e~*^). This expectation was introduced by Bui and Keating in |BuKej (following |GHKj ) as the 
random matrix counterpart of a key term contributing to the mean values of certain Dirichlet L- 
functions in the Katz-Sarnak orthogonal family. The issue is the large n and large X behavior of 
this average. The question can be resolved with the help of (jl75p . One observes that 



(184) 



E 



50(2n) 




^ det {(t>j^k + 4>j+k)lkLo 



and that the symbol (|182p is of FH type with a single a-singularity at zq = 1 and ao = k. Directly 
from (the precise version of) the asymptotic formula ()175p (see Theorem 1.25 in [DIKlj ). one 
obtains, for X large, as n — )• oo 



(185) 



E 



SO{2n) 




Gil + k) 



r(l + 2A;) 



G{l + 2k)r{l + k) J \e-'E log X 



2n 



fc(fc-i) 

2 



where 'Je is Euler's constant. Formula (jl85p is precisely the asymptotic form conjectured by Bui 
and Keating in |BuKe| . 



10 Block Toeplitz matrices 

We now consider block Toeplitz matrices Dn{f) of the form Dn{^) = det((/3j_fc)"^^Q where ifq € 
M(r, C) are the Fourier coefficients of an r x r-matrix- valued integrable function y5(e'^) on . The 
earliest asymptotic result for such determinants is the analog of Szego's Theorem 1 due to Gyires 
in 1956 [Gyj , who showed that if ^p{z) is continuous and positive definite for \z\ = 1, then 

(186) lim i log L>„ ((/?) = — / logdet (^fe*^) 

Block Toeplitz matrices arose in the Ising model in the following way. In their famous paper 
[LeYaj on the Lee- Yang Theorem, asserting that the zeros of the partition function Zp^^h (cf. Re- 
mark 14.11 in Section U]), for the Ising model in the presence of a magnetic field, all lie on the 
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imaginary /i-axis, Lee and Yang also proposed, but did not prove, expressions for the free energy 
Fh = lini|A|^oo ~kB T log Z/i^h and the magnetization M/j = lim„_^oo (^i,! ai^i+n) 2 at the particular 
imaginary value of h, = ^ ztt. In terms of the variable p = q-'^^/^bT ^ ^j^g zeros of ^a,/i lis on 
the circle \p\ = 1. As noted in [LeYaj . is precisely the density of the zeros of Z^^^ on {\p\ = 1} 
at p = —1, in the limit as |A| — t- 00. In [McWu2| . McCoy and Wu present a derivation of these 
expressions for and M^, h/k^ T = ^iir. Regarding M^, for this value of h, they first show how 
to express the correlation function (fxi^i ai^i^n) in terms of a 2 x 2 block Toeplitz determinant 
Then they show by certain ingenious manipulations how to evaluate Bn in terms of a determinant 
which is a product of two scalar Toeplitz determinants to which SSLT applies. In this way they 
are able to compute Mh explicitly and recover the formula in |LeYa| . This is the first example of 
an SSLT-type result for block Toeplitz matrices. 

The general theory of SSLT for block Toeplitz matrices begins with Widom's paper [WidSj . 
Under certain smoothness assumptions on the r x r-matrix- valued function 99(6*^), Widom shows 
that 

(187) Y{ip) = lim D,((^)/exp logdet ^(e'^) do] 

n->oo IZTT Jq V / J 

exists, and in [Wid3] he shows that 

(188) Yi^) = det[T{^)T{^~')] 

where T{ip), T((/3~^) are the Toeplitz operators associated with ip and (p>~^ respectively (cf. ()76p ). 
Part of the proof of (|188p consists in showing that T{(p) T{ip^^) — 1 is trace-class in so that the 
RHS is well-defined. In the scalar mentioned above, Widom uses the Helton-Howe-Pincus 

formula to show that (fTHH]) reduces to the standard result Y{ip) = eS^=i ''(^"S'p)k (logv')-*. As noted 
above (see (f73]) et seq), in |Wid4] . Widom introduced a "standard" pathway to SSLT. The method 
extends with only a change of notation to the block Toeplitz case, and so Widom obtains a new 
direct proof of ()187p . (jl88p with no more effort than in the scalar case. The same is true for the 
method of Basor and Helton in [BasHelt] . 

Block Toeplitz matrices occur in many physical situations, for example, in the study of the Ising 
model with next-nearest-neighbor interactions, and also, more recently, in the theory of the dimer 
model [BasEhrj and in the study of the entanglement spectrum in quantum systems (see [VLRK] . 
[ JK] ) . The use of the general SSLT for block Toeplitz matrices in physical applications poses new 
challenges vis a vis the scalar case. The new difficulties are due to the fact that formula (jl88p . as 
elegant as it is, is hard to evaluate in concrete calculations in the case of a matrix symbol (/?. The 
determinant on the right hand side of (jl88p is the Fredholm determinant of an infinite matrix, and 
even for 2x2 matrix functions if an effective evaluation of Y(ip) is a highly nontrivial enterprise. 
Indeed, the authors are aware of only three situations when such effective evaluations have been 
made so far. The first situation arises when, as in the work of McCoy and Wu mentioned above, 
it is possible, by using an ad hoc technique, to express the block Toeplitz determinant in question 
in terms of scalar Toeplitz determinants. The work of Tanaka, Morita, and Hiroike [TMH] . and of 
Bottcher [Bott2| on the Ising model with the next-nearest-neighbor interactions provides another 
example of this sort. In this example, the relevant matrix symbol <p is described by the formula 

(189) ip{z) = {I-zR)a{z){l-z-^S) 
where a{z) is a diagonal matrix function, 

'X{z) 
X{z-^), 



(190) a{z) 
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with X{z) being an elementary algebraic function, and the constant matrices R and S are expressed 
as follows: 



^-Ko '^"{o :)^-' 



In these expressions, U is an invertible 2x2 matrix, and q, /3 are complex numbers, satisfying the 
conditions: \a\ < 1, \/3\ < 1. Note that these inequalities guarantee that the first matrix factor in 
()189p is holomorphic and invertible inside the unit disc while the last factor is holomorphic and 
invertible outside of the unit disc. The function X(z) (an analog of Onsager's symbol ()25p ). the 
matrix U, and the numbers a and /3 are determined by the physical parameters of the problem, the 
temperature T and the anisotropy coefficients. As in the usual Ising model, there exists a critical 
temperature Tc, and the behavior of the determinant Dn{f) depends on whether T > T^ox T < Tc- 
Bottcher shows that in the case T > Tc (the paramagnetic phase) the Wiener-Hopf factorization 
of the matrix symbol ip~^{z) has nontrivial partial indices which implies that the operator T{ip~^) 
is not invertible. In view of (jl88p . this implies that the factor Y{ip) is zero. Moreover, it is easily 

p27r . 



seen that J^^ logdet ip{^e jdO = 0, and hence for T > 

lim Dn{v) = 



which has the physical interpretation that there is no spontaneous magnetization in the paramag- 
netic phase. In the ferromagnetic phase, that is when T < Tc, an efficient direct analysis of Widom's 
factor Y{(p) is not apparent, and Bottcher develops an ad-hoc approach which uses the specifics 
of the symbol (11890 . Indeed, in his preprint, by means of ingenious operator-algebra arguments, 
Bottcher obtains the formula 



(191) 



lim 



-Dn-l(y) 

DJa) 



{deiU)-\l-a(3)- 



U11U22 



U12U21 



A+(/3)A_(l/a) A+(a)A_(l//?) 



where A = A+A- is a Wiener-Hopf factorization of A (which, as Bottcher shows, is possible in the 
case T < Tc). Hence the evaluation of Dn-i{y^) as n — >• 00 is reduced to the asymptotic evaluation 
of the scalar Toeplitz determinant Dn{X), for which the standard (scalar) SSLT applies. In the 
critical case, T = T^, the function A(z) acquires a Fisher-Hartwig type singularity. Bottcher then 
provides a modification of his previous analysis for this case and, using in addition various results 
obtained earlier with Silbermann (see Section [6]), he then derives the leading asymptotics of Dn{^) 
in the critical case as well. 

The two other types of block Toeplitz determinants for which the Widom pre-factor Y(ip) can 
be evaluated are not reducible to the scalar case. The first one corresponds to the matrix symbols 
(p with a Fourier series that is truncated on at least one side. This class was singled out by Widom 
himself in [WidS] . For such matrices if, the evaluation of Y{ip) is reduced to the evaluation of a 
finite dimensional determinant (in fact, a block Toeplitz determinant of a finite fixed size). This 
result of Widom has been used in the recent paper |BasEhr| of Basor and Ehrhardt devoted to 
the dimer model. The matrix function (p appearing in the dimer model is not originally of the 
truncated form, but it is truncated up to a scalar algebraic factor. Basor and Ehrhardt show that 
this factor can be accounted for by means of skillful algebraic manipulations which transform the 
original quantity Y{ip) to one with (p replaced by its truncated part. 

The second type of block Toeplitz determinants that are not reducible to the scalar case, but for 
which the large n limit is still computable, are determinants with algebraic matrix symbols. Such 
determinants appear in the theory of quantum entanglement when one evaluates the von Neumann 
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entropy, which is a fundamental measure of the entanglement of a subsystem with the rest of the 
quantum system. For instance, as shown in |VLRK| . [JKj . the evaluation of the von Neumann 
entropy of a subsystem of neighboring spins in the XY quantum spin chain is equivalent to the 
asymptotic evaluation of the block Toepltz determinant Dn{'-p) whose matrix symbol ip is given by 
the formula (for more details and references see [JKj ) . 



Here, zi ^ Z2 are complex nonzero numbers not lying on the unit circle and determined by the 
physical data of the model (the anisotropy parameter and the magnetic field) and // is a free complex 
parameter lying outside of the interval [—1,1]. The asymptotic evaluation of Dn{(p) is performed in 
[IJKj and it is based on yet another formula of Widom which he obtained in [WidSj . This formula 
is concerned with block Toeplitz determinants depending on an external parameter, say fi, and 
is used to evaluate the leading asymptotics of the ^-logarithmic derivative of the determinant in 
terms of the Wiener-Hopf factorizations of the matrix function (p~^ , i.e. in terms of the it - matrix 
functions u± and v± defined by the equations 

(193) (f^^{z) = U-^{z)u-{z) = V-{z)v-^{z) u_(oo) = w_(cxd) = /, 

where, as usual, /+ (/_) means that the matrix function is invertible and analytic inside (outside) 
of the unit circle, respectively. Under some natural smoothness assumptions (see |Wid5] . Theorem 
4.1), Widom's formula reads: 



(In the case of symbols analytic in an annulus including the unit circle, formula (|194p was re-derived 
with a quantitative estimate of the error term in |IJK] via the Riemann-Hilbert scheme (see also 
[IMMj l.l Widom used this asymptotic formula in the proof of his main result, vis., the asymptotic 
relation ()187p . A curious fact is that, up until the work [IJKj, equation ()194p . apparently, had never 
been used directly for the asymptotic evaluation of the determinants of block Toeplitz matrices. 
And for a good reason: the efficiency of formula (I194p relies on the effectiveness of the Wiener- 
Hopf factorizations (jl93p of a given matrix valued function (p{z). Usually, this is a transcendental 
problem equivalent to the solution of a system of singular integral equations on the unit circle which 
can seldom be effected in terms of known special functions. However, as was observed in [IJKj . in 
the case of algebraic matrix functions (/'(•z) one can take advantage of the algebro-geometric method 
that had been developed in the late 70s to the early 90s in soliton theory (see e.g. [BBEIM] ). Using 
the Riemann-Hilbert version of this method exploited in [DiZ] (see also |DKMVZ1] . |DKMVZ 2]). 
the authors of [UK] were able to construct the Wiener-Hopf factorizations of the symbol (I192p in 
terms of Jacobi theta- functions. This in turn, with the help of Widom's formula (jl94p . yielded the 
following asymptotics of the corresponding Toeplitz determinant: 



(192) 




and 




(194) 




(195) 




n — oo. 
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Here 



oo 



fe=— oo 



is the Jacobi theta- function, cr = ±1, and r is the modulus of the ehiptic curve 



(196) w\z) = {z- zi){z - Z2){Z - z^^){z - z^^). 



The concrete choice of a depends on the way the branch points z^ ^ ^"^^ located with respect to the 
unit circle (see |IJK] for details). We refer the reader to the survey [IK] for more on these results 
and on their applications in the theory of quantum entanglement. The asymptotic formula ()195p 
was extended in [IMMj to the case of more general integrable spin chains introduced in [KM] . The 
relevant symbol has the same matrix structure as p92|) . but with scalar function "il^iz) now given 



where ^{z) is a polynomial of degree 2m. The analog of the formula (jl95p in the case m > 1 involves 
a hyperelliptic curve instead of an elliptic curve, and instead of the Jacobi theta-function one needs 
2m — 1 dimensional Riemann theta-functions. Independently of [IJKj and [IMMj . similar results 
for block Toeplitz determinants with algebraic symbols appearing in the theory of the Gelfand- 
Dickey hierarchy, were obtained, again with the help of the algebro-geometric approach, in [Caf]. 
It is worth noticing that the method of [IJKj can be applied to the problems considered above by 
Bottcher and Basor-Ehrhardt. Indeed, the symbol of the determinant appearing in the work of 
Basor and Ehrhardt is in fact algebraic and hence can be factorized within the algebro-geometric 
approach, while the symbol of Bottcher's determinant is already in the form uj^U-. The reverse 
factorization, i.e. representation of the symbol ()189p as V-v^^ involves an elementary algebraic 
operation. We note finally that Widom's formula (|194p can be used to extend Bottcher's results 
to symbols of the form ^{z^ = R{z)a{z)S{z~^), where a{z) is a diagonal matrix function and R{z) 
and S{z) are polynomial matrix functions invertible for all z inside the unit circle. 

11 Double-scaling limits 

Of all the challenges in Toeplitz theory that arise from the Ising model, the deepest is the issue of 
the double-scaling limit (or in physicists' parlance, simply the scaling limit), T — )• Tc and n — )• oo. 
We have already mentioned the importance of this problem in physics. At the purely mathematical 
level, the issue is the following: Suppose one has a Toeplitz determinant Dnif) with a symbol that 
depends on some external parameter, say t, 93(e*^) = ^pt{e'^). For t > 0, (pt is regular, but at t = 0, 
(ft has a (Fisher-Hartwig) singularity, and hence, as we know, the nature of the asymptotics of 
Dn{ipt) as n — )• oo is different for i > and t = 0. This raises the general question: How does one 
describe the transition from the one kind of asymptotics to the other as n — )• oo and t ^0? 

In [WMTB] ■ Wu, McCoy, Tracy, and Barouch discovered a remarkable scaling relation for the 
Ising model. They defined the scaling functions 



by 



(197) 




(198) 



G±{r) 



lim 

£,m— >oo,t— >=fO 



1 - 
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in the double-scahng hmit 



(199) 



sinh (2 Ji/kB T^) f + sinh (2 Ja/A;^ 
sinh (2 Ji/kB T^) + sinh (2 Js/Zcb Tc) 



|f| = r, fixed 



Here it refer to T > Tc and T < Tc respectively. They then showed that G±{r) could be expressed 
in terms of a solution of the Painleve III equation 



(200) 



as 



(201) 



G±(r) 



cPt] 
1¥ 



l=Fr? 



2r/(i) 



1/2 



drj 

de 



exp 



1 dr] o 

+ V 

Ode ' 



V 



ji-v')'- iv') 



r/2 



/\2 



dO 



with the boundary condition r]{6) ~ 1 — |: Xo(2^) as ^ — t- oo, where Kq denotes the modified Bessel 
function. The above calculation was the first derivation of an explicit scaling law for a two-spin 
correlation function for any model in statistical physics. As noted before, the scaling functions 
G±{r) are believed to be universal for a wide class of two-dimensional models with short range 
interactions. 

A crucial issue for |WMTB] scaling theory was to show that formula (|198|) matches the critical 
1/4 behavior (jl06p as r — )• 0. In order to do this, the authors needed to know the oo -H- 
connection formulae for the Painleve functions r]{0). In |WMTB] they extracted this information 
from the unpublished thesis of Myers [My] (see below). In a later paper in 1977, McCoy, Tracy, and 
Wu |MTWj derived oo o connection formulae for a two-parameter class of bounded solutions 
r]{0;v, A) of the general Painleve III equation 



dFr] 1 f dr]\ 1 dr] 2v 



r]\dej Ode e 



V 



We now describe in more detail what was done in [MTWj . restricting our discussion to the case 
u = which is relevant to [WMTBj . 

Consider the one parameter family ri{9) = r]{9;0,X) of the solutions of equation ()200p defined 
by the asymptotic condition 



(202) r]{0) ^l-2XKoi2e), 6 

In [MTW] . it was shown that, for < A < l/vr, as r — )• 0, 

1 



oo. 



(203) 



r/(r/2) = Br"" 1 



16 



-B-^(l 



(7 



2^2-2<7 



+ 0(r2 



where the constants B and a are functions of A given by the explicit formulae: 



(204) 



a = (t(A) = — arcsin(7rA), B = B{a) = 2 



IT 



3, r((l-a)/2) 
r((l + a)/2)- 



In |My| , Myers derived the asymptotics for r/(r/2; 0, l/vr) as r | 0, 

1' 



(205) 



rj 



r 

-;0, - 

2 TT 



log -+1E 
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where 7^; is Euler's constant. In [ MTWj the authors recovered (|205p by analyzing formally the 
limit A — )• I/tt in ()203p . Proofs of (|205p were given later in [ WidSj and, using the Riemann-Hilbert 
method, in [Nilesj . In turn, the asymptotic formula (j205p yields the estimate: 

(206) G± ~ const r"^/^ 

In order to see that this estimate matches equation ()106p . it is enough to notice that (|198p and 
(fT99]) imply that 

(207) ~ 2'/^Gi{r)\t\'/' = 2'/^Gi{r)r'/\-'/\ 

The main term in equation ()106p . up to the multiplicative constant, follows now directly from ()206p . 
The evaluation of the constant in (j206p that matches the transcendental constant A in (jl06p was 
done by Tracy [T] in 1991. 

The Painleve functions, PI, . . . , PVI, (see [Ej) were discovered and analyzed intensively at the 
beginning of the 20*^^ century until the First World War, and then fell into a period of latency. They 
reappeared in theoretical physics quite unexpectedly in 1965 in the work of Myers |My[ mentioned 
above, who showed that the scattering of electromagnetic radiation from a strip in the plane could 
be expressed in terms of a solution of the Painleve III equation. The paper [WMTB] was a seminal 
event in mathematical physics, followed soon after by the landmark paper of Jimbo et al [JMMS] 
on the impenetrable Bose gas, and the Painleve equations have emerged over the years as the 
core of modern special function theory, with applications across the board in mathematics, physics 
and engineering (see, e.g., |FIKN] ). In PDE's, in particular, the Painleve equations (specifically, 
Painleve II) appeared for the first time in the analysis of Ablowitz and Segur |AS| of self-similar 
solutions of the Korteweg-de Vries equation. This pioneering development opened the door to 
understanding the crucial role that the Painleve equations play in the theory of integrable dynamical 
systems (see |AC| . [FIKN] ) . The work [MTW] demonstrated for the first time the possibility to 
derive explicit connection formulae for families of solutions to a Painleve equation, a property which 
had been previously known only for the linear ODEs of hypergeometric type. Similar connection 
formulae were obtained for the second Painleve equation in {AS] (only formulae for the amplitude) 
and in |SA| (complete formulae for the amplitude and the phase). The papers |MTW| . |AS| . and 
|SA| generated a surge of activity in the global asymptotic analysis of the Painleve equations (see 
again |FIKN] for more history and the state of the art in the area). We view these developments as 
one more example of the way in which questions arising in the analysis of the Ising model gave rise 
to developments in Toeplitz theory, and in this case, through Toeplitz theory, to the global theory 
of the Painleve equations. 

Recently in [CIK], the authors considered, in particular, the double-scaling limit for the diagonal 
correlation function {ai^iai^n,i+n) = Dn{f ding) (see (f3T]l . (|32]) ) for the Ising model as n — )• 00 and 
T t Tc. Since this correlation function is given by a Toeplitz determinant, it is possible to consider 
the double-scaling limit rigorously using Riemann-Hilbert/steepest descent methods. The authors 
define the scaling parameter t in the same way as in (jl98p . i.e. e~* = fcons, so that 1 10 corresponds 
to T t As noted above (/^diag is regular for t > and has a Fisher-Hartwig singularity for t = 
at zq = 1, with ao = and /3o = ¥'diag(6') = e*^''-^)/^ (see (HM^i ). The main result in [CIK] is 
the following: There exists a finite (perhaps empty) set A C (0, 00) and a (small) number to > 0, 
such that for n — t- 00 and for all < t < log -Dn(v'diag) has a uniform expansion of the form (see 



r /, r \ 
1±. ln-+7i. 
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[UIK] Theorem 1.4) 

1 . .,,1 Z-^"* cr(j;) + 1/4 

-dx 



1 1 /" 

log (cri,i ai+n,i+n) = log -Dn(v^diag) =7 log (l - e~^*) - - log 2nt + I 

(208) ^ . ^ / 

^1\ ^ /3 



provided dist (2nt, A) > J > for some 5 > 0, and the path of integration in (j208p does not 
intersect with A. Here cr(a;) is a solution of the Painleve V equation 



satisfying the boundary conditions 
(210) a{x) = 



-l/4 + 0(3;lnx), x 0, 



2ix 



X 



-ie-(l + 0(i)), x^+oo. 



In addition, the following integral identity for the function ct(x) holds true (see (1-31) in [CIK]): 

^d.-i,og.„ = -logG(i)o(f' 



,211) r^M±izi..+ r^M<,.-i,og.„=-iogGfiVfli. -o^A 



which together with ()208p and the definition ()198p of the scaling functions G±{r) yields an alter- 
native to the representation (j20ip of G-{r) in terms of the fifth Painleve transcendent a{x): 



(212) G_(r)=exp 



2r 



This formula was obtained previously in 1980-81 by Jimbo and Miwa [ JMj . The equivalence of 
(j212p to the original PHI- formula (j20ip was established in [MPJ . Hence the uniform asymptotics 
(j208p includes the double-scaling limit 2nt = 2r — )• const (away from A), n — )• oo, in |WMTB] as a 
special case. 

Remark 11.1. In |CIKj . the authors derive (|208p for symbols which are more general than (/^^iag, 
specifically for symbols that have a general FH singularity at f = with Kao > —1/2, /?o £ C, 
ao ± /3o / -1, -2, . . . (see (1.7) in [CE]). 

Remark 11.2. In |JMj . the authors used a remarkable relation (which they themselves discovered) 
between the diagonal Ising correlations and isomonodromy deformations of a certain 2x2 Fuchsian 
system. This discovery laid the foundation for the Riemann-Hilbert method in the theory of 
correlation functions, random matrices, and Toeplitz and Hankel determinants (see |DIZ] for more 
on the history of these developments). Using this link to Fuchsian systems, it was shown in |JM] 
that the diagonal correlation function {(Ji^i o"i+„^i+„), as a function of the variable 

s = e~^* = ( sinh - — I- sinh - — ^ | 
V kBT kBTj 

is expressed in terms of a solution to the Painleve VI equation. The derivation of ()212p in |JM] 
is done by formally performing a scaling transformation taking the inverse monodromy problem 
associated with this Painleve VI to the inverse monodromy (Riemann-Hilbert) problem associated 
with the fifth Painleve function a{x). In the rigorous setting of [CIK] . the latter problem appears 
as a parametrix during the implementation of the nonlinear steepest-descent method. 
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Remark 11.3. The Painleve VI description of the diagonal correlations in [JM] was the first example 
of the appearance of the Painleve functions in the theory of Toeplitz determinants before the large 
n limit is taken. One realizes now that nonlinear differential equations of Painleve type are always 
present in situations where the symbol (/'(•z) satisfies the condition: 

(213) ^/^^ ^ = rational function. 

dz 

This fact is quite easy to see within the Riemann-Hilbert formalism (see, e.g., discussion in [IT W] ) . 
These "finite n" Painleve representations are not universal: the type of the nonlinear ODE depends 
strongly on the structure of the rational right hand side of ()213p . i.e. on the symbol. 

Double-scaling problems for Toeplitz determinants occur in many different areas. For example, 
in combinatorics, let vr = (7r(l), 7r(2), . . . ,7r(A^)) be a permutation of the numbers 1,2, ... ,A^. If 
1 < ii < ^2 < • • • < ^fc < and 7r(ii) < 7^(12) < • • • < vr(zfc), we say that '/r(ii), . . . ,-K{ik) is an 
increasing subsequence in vr of length k. Let i]\f{iT) denote the maximum length of all increasing 
subsequences in vr. With uniform distribution on the permutations, one is interested (see [BDJj ) 
in the distribution function 

PN{n) = Prob {£Ar(7r) < nj = — 

as A^, n — )• 00. It turns out [Ges] that the Poissonized version </)n(A) of the process PN{n) 

(214) <i>n{\) = ^"^]v! ^ ^ ° 

can be expressed in terms of a Toeplitz determinant 

(215) <PnW = e-^ Dni^x) 

where (l)\{9) = g^^™*^^. In order to recover pN{n) as A^, n — t- 00 from (|214p (|215p one must 
analyze the double-scaling limit for the determinant as A, n — t- 00. The main result in [BDJj is the 
following: 



(216) ^,i^„^P„b (^£^<tj=f,«) 

where F2{t) is the Tracy- Widom distribution |TWj for the largest eigenvalue of a random matrix 
in the Gaussian Unitary Ensemble. The distribution F2{t) can be expressed in terms of Painleve 
functions as follows: 

(217) F2(t) = e~ii°° 
where u{s) is the unique solution of PII, 

(218) u"{s) = 2u^{s) + su{s) 
with 

m(s)~— Ai(s) as s — +CXD. 



52 



Here Ai (s) is the standard Airy function. On the other hand, F2{t) can be expressed as a Predholm 
determinant 



(219) F2{t) = det{I - K^^^^) 

where -f^'Xiry trace-class operator with kernel 



^is) , , Ai(x)|Ai(y)-Ai(y)AAi(x) 



X -y 

acting on L^(s,oo). In random matrix theory, this determinant is the probability that the interval 
(s,oo) contains no eigenvalues in the edge scaling limit in the Gaussian Unitary Ensemble. 

The result in [BDJj was the first of many results linking problems in pure and applied math- 
ematics and in physics to the Tracy- Widom distribution F2{t) and its symplectic and orthogonal 
analogs (see, for example [Dei2| ). 

As already noted on several occasions, in his work on impenetrable bosons [Lenlj, Lenard was 
led to consider Toeplitz determinants with symbols of type |e*^ — e'^^ | |e*^ — e*^^| (see (|125p above) 
that vanished on S^. In |Len2j he then considered symbols with any finite number of such a- type 
FH singularities. Such symbols vanish at only a finite number of points, and in [Wid6j Widom 
began considering symbols which vanished on a full interval. Let (/9^(e*^) denote the characteristic 
function of the interval (//, 27r — //), < /x < vr. In a virtuoso calculation Widom showed that, for 
jj, fixed, as n — )• 00 

(220) log Dni^^) = r? log cos | - ^ log (n sin |) + co + o{\) 

where cq = ^log2 + 3C'(— 1) and C,(^z) is the Riemann zeta-function. A few years later in 1976, 
Dyson |Dy4| returned to the problem (cf. discussion following (j59p above) of computing the prob- 
ability Pg that there are no eigenvalues for a random matrix in the interval (0,2s/7r), in the bulk 
scaling limit for the Gaussian Unitary Ensemble. He showed that as s — )• 00, logP^ has a full 
asymptotic expansion 

(221) logP, = -^-^logs + ao + — + ^ + ... , 

2 4 s 

Dyson identified all the constants ao, 01,02, • • • , and it turns out that oq, which is of particular 
interest, is precisely the constant cq arising in Widom's expansion (j220p above. Dyson arrived at 
the identification oq = cq by noting first that the probability Ps is given by (see, e.g., |Meh] . [DeiS] ) 

(222) Ps = det(/ - Ks) 
where Ks is the trace-class operator with kernel 

(223) i^.(x,,) = !i^^ 

7r(x - y) 

acting on L^(— s,s). A simple calculation shows that 

(224) D^{^,) = det (5,, - ^^^^(^•-^) '\ 

V 7r(j-fc) Jo<j,k<n-l 
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and so for fixed s > 0, 

(225) lim Dn(ip2s] = Ps 

and hence, if the error term o(l) in ()220p was uniform in the double scaling limit, n — )• oo, ^ — )• such 
that fj,n = 2s, one could conclude from (|220p (j225p that ao is indeed equal to cq. The uniformity of 
the error term o(l), however, remained open. There are now three proofs that ao = cq. The first two 
were given simultaneously and independently by Ehrhardt and Krasovsky, and the third was given 
a little later in [ DIKZj . The proofs by Krasovsky and by Deift, Its, Krasovsky, and Zhou proceed 
by following Dyson's observation (j225p and verifying that the o(l) error term is indeed uniform. We 
refer the reader to |DIKZ] for some history of the problem P^, s — )• oo, and, in particular, for the 
relevant references to the work of Krasovsky and Ehrhardt. A similar asymptotic problem arises 
when instead of the sine kernel determinant (1222p one considers the general so-called confluent 
hypergeometric kernel determinant which depends on two extra parameters (and is related to the 
so-called Bessel kernel determinant and also to the determinant (j222p for particular choices of the 
parameters); another related problem is the asymptotics for the Tracy- Widom distribution (j219p 
as s — >• — oo. These problems can also be represented as double scaling problems for Toeplitz (and 
also Hankel) determinants, and such representations were crucial for finding complete solutions for 
the asymptotic problems at hand. We refer the reader to |Kr| for a discussion of the results and 
for the references. 

Dyson's identification of 01,02,... in |Dy4| involves an ingenious application of inverse scat- 
tering theory for Schrodinger operators on the half- line < x < 00. He begins by noting that 
Ps = det(/ — Ks) can be factored 

(226) det(/ - Ks) = D+{s)D^{s) 

where D±{s) = det(I — /±)l2(o,s) are Toeplitz+Hankel type determinants on L'^{0,s) with kernels 

(227) f±{x,y) = - 

vr 

He then considers the functions 

(228) W±{s) = -2^ log D±{s)-1 

as potentials for Schrodinger operators H± = —-^ + W±{s) on L^(0,oo) with appropriate Robin 
boundary conditions at s = 0. Using the Gelfand-Levitan formalism (see, e.g., [Fad] ) j Dyson is able 
to identify the spectral measures p±{X)dX for the operators H±. He then obtains the scattering 
phase functions e*^^. In the Marchenko formalism (see, e.g., [Fad] ) one utilizes the functions e*^^ 
to express the potentials W±{s) in terms of Fredholm determinants of the form 

(229) w±{s) = -^(^s±^^ _2^1ogA±(s), s>^ 

where A-i-(s) = det{I — F±)i2(^s,oo) and F±{x,y) are certain explicit kernels which decay as s — )• 00. 
Equating ()228p and ()229|) . and integrating twice, one obtains the following identity for s > ^: 

(230) logP, = log(/ - Ks) = log D+D^ = - i>,2 _ ^ j^g (^2 _ ^ V + i^g (A+(s)A_(s)) 



sin(x ~ y) _^ sin(x + y) 
X — y X + y 
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(the terms linear in s drop out). As F±{x,y) decay as x,y oo, it is clear that (|230p gives rise 
to a full expansion of the form (|22ip . where the coefficients can be computed in terms of the trace 
powers trF^. 

Dyson is able to compute e^^^ , and hence F±, explicitly, because of the following fundamental 
fact from inverse scattering theory. The spectral weights have the form p±{X) = \/A/ (7r|v7±(\/A)p), 
where ^±{k), k = ^/X, are certain functions constructed from the Jost solutions J±{s,k) for H±, 
H±J±{s, k) = k'^J±{s, k), J±{s, k) ~ e^^'^ as s — )• oo. The functions 'p±{k) are analytic in {Qfc > 0} 
with prescribed asymptotics ask ^ oo, and in the cases at hand, have no zeros in > 0}. Hence, 
by standard arguments, ^p±{k) are determined by their absolute value \^{k)\ for "^k = 0. Thus p±{X) 
determine ^±{k). However, e^^^ can be expressed in terms of as e^"^^^^^ = i'^±{k) /\^±{k)\^ and 
hence p±{X) determine the scattering phase functions, and hence, eventually, F±. 

These considerations bring to mind the Borodin-Okounkov-Geronimo-Case formula (j43p in 

which the Szego function ^^{z) = exp f^^^f log 99(6*^) dO^, which is analytic in \z\ < 1, 

plays the role of the functions ip±{k). Indeed, |P(e*^)p = ip{e^^), which is the spectral measure 
for the associated Cantero-Moral- Velasquez (CMV) unitary operator acting on £^ (see [Sim2]). 
Thus we may think of the Toeplitz determinant D^i'-p) as an expression in the "Gelfand-Levitan 
formalism" . On the other hand, we may think of the Predholm determinant on the RHS of (j43p as 
an expression in the "Marchenko formalism". Indeed, from ()44p . b = T) /iPp, which we may view 
as a scattering phase function b = e^*^ for the CMV operator, and c is just b = e"^*''. These con- 
siderations are very suggestive and invite further investigation. One notes that the original proof 
of (P3|) by Geronimo and Case in [GerCase] was developed in the context of scattering theory. 

Another example of a double-scaling limit is the following. Recall that the density matrix for 
the impenetrable bosons is given by pn,l{0 = z^Ni'^'^i/ L), where Rn is the Toeplitz determinant 
(jlSOP with symbol (I13ip . For definiteness, fix the length scale so that L = N. As we discussed 
in Remark 16.11 in Section [6l of particular interest is the double-scaling limit N — )• 00, ^/N — )• 0. 
In [Lenlj Lenard considered this limit p{^) = limAr_s.oo j/RnC^t^^./N) with ^ fixed. He obtained 
an expression for p($^) in terms of a Fredholm determinant with explicit kernel related to the sine 
kernel (1223p . He also obtained another representation for p(^) in terms of a Fredholm minor for 
the kernel 2 "°"I^I^V^ . 

TT{x-y) 

In order to complete his analysis of the limiting momentum distribution in the gas, Lenard 
needed to know the large ^ behavior of p(^), but he did not derive this behavior in [Lenlj . The 
derivation of this limiting behavior was carried out by Vaidya and Tracy [VaiTrj who found that 

(231) piO = (1 + ^ (-^ 2vre - ^) + • • •) , p^ = vreV22- V3^-6, 

as ^ — )• 00, where A is again Glaisher's constant (jl07p . Note that if one evaluates Pn,n{0 = 
■^RNi^TT^/N), N,^ = ^ —f- 00 using ()137p ()138p . and assumes that the formulae remain valid as 
t = 2iT^/N — )• 0, then one obtains the leading term in (|23ip . 

The double-scaling limit — )• 00, ^/N — )• corresponds to the Toeplitz determinant with two 
FH singularities merging (at z = 1) with ai = a2 = |, /3i = /32 = 0. A general analysis of 
the double-scaling asymptotics of a Toeplitz determinant with two merging FH singularities (in 
particular, of Lenard's symbol (|130p ) is given in [CK j. 

Dyson also considered [Dy5| the following scaling problem associated with Toeplitz determi- 
nants. Let 

(232) /^{z,t) = dQl{I - zkt) 
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where < z < 1 and Kt acts on L'^{—t, t) with kernel Kt{x, y) = — ^^^fr^y^- With the extra vr in the 
sine function (cf. Kg in ([223])), A(l,t) is now the probabihty that there are no eigenvalues in the 
interval (0,2t). Of interest here is the double scaling limit for A(z,t) as t — >■ oo and z t 1- Dyson 
analyses the problem by interpreting /S.{z,t) in terms of a Coulomb gas at inverse temperature 
/3 = 2. In this interpretation, z corresponds to an external potential v via the formula 1 — z = e~^'". 
Then 

where 

(234) Z;3 = J]e-^(^+^") 

is the partition function of the gas with external potential v applied to charges in a fixed interval 
of length 2t. The sum in (j234p is to be interpreted as an infinite-dimensional integral over all 
configurations of the gas, W represents the sum of all Coulomb interactions Wjk = — log \xj — Xk\, 
and for any configuration, n is the number of charges in a fixed interval of length 2t. Dyson does 
not define the sum rigorously, but he uses it as a heuristic guide in his calculations. The main result 
in the paper is the calculation of oscillatory factors in A(z,t) as t — )• oo: These factors turn out to 
be so-called "genuinely non- linear oscillations" and are described by Jacobi elliptic functions. For 
< z < 1, < A < 1, the symbol 

(235) fz,\{e'^) = X[7rA,7r(2-A)] + (1 " 2)X[-7rA,7rA] 

gives rise to the Toeplitz determinant 

(236) DM.,x) = det (s,, - z'^^^^^t_3\ 

V <3-k) Jo<j,k<n- 
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and so for fixed z, t 
(237) llm Dn(f,2t) =Aiz,t). 



Thus Dyson's double scaling limit z f 1, t — )• oo can be viewed as a triple scaling limit for the 
Toeplitz determinant Dn{fz,2t/n) with z'l l,t ^ oo, and n — >• oo. Note finally that fz^\ is precisely 
the kind of symbol that arises in the Toeplitz eigenvalue problem (|178p . 

Lenard ends his 1972 paper |Len2] with a hope and a prophecy: "It is the author's hope that a 
rigorous analysis will someday carry the results to the point where the true role of the zeros of the 
generating function will be understood. When that day comes a capstone will have been put on a 
beautiful edifice to whose construction many contributed and whose foundations lie in the studies 
of Gabor Szego half a century ago" . 

Almost 100 years have passed since Szego published his first paper on the asymptotics of Toeplitz 
determinants. We certainly know more now than in 1972, but many new problems continue to arise, 
and it is still too early, alas, to set the capstone. 
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